• Economic and Environmental Geology
• /
• v.22 no.3
• /
• pp.267-276
• /
• 1989

This paper presents results of interpretaton of gravity data by iterative nonlinear inversion methods. The gravity data are obtained by a theoretical formula for two-dimensional 2-layer structure. Depths to the basement of the structure are determined from the gravity data by four interative inversion methods. The four inversion methods used here are the Gradient, Gauss-Newton, Newton-Raphson, and Full Newton methods. Inversions are performed by using different initial guesses of depth for the over-determined, even-determined, and under-determined cases. This study shows that the depth can be determined well by all of the methods and most efficiently by the Newton-Raphson method.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.379-382
• /
• 2003

In gravity data correction process, mass effect of the upper part of base level is removed with Bouguer density. Usually, Bouguer density is estimated as a mean density in the field area. But, this may causes a serious problem when ore body is in the area. To overcome this problem, we tried to apply a new method mixing up mass corrections and inversion (3DGTI). 3-D Gravity Terrain Inversion (3DGTI) includes information of topography and distribution of Bouguer density. For this method does not remove the mass effect above base level, it is no longer useless to use Bouguer density. Numerical model tests have shown that the 3DGIT successfully retrieves the anomalous subsurface density distribution of both surface and deeper layers. Model tests shows that this method shows better results than those of conventional one, especially when main target is ore body. The inversion result well delineates the three-dimensional shape of the intruded granite body and basement.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.479-489
• /
• 1994

In this paper, the iterative least-squares inversion method is used to determine shapes and density contrasts of 2-D structures from the gravity data. The 2-D structures are represented by their cross-sections of N-sided polygons with density contrasts which are constant or varying with depth. Gravity data are calculated by theoretical formulas for the above structure models. The data are considered as observed ones and used for inversions. The inversions are performed by the following processes: I) polygon's vertices and density contrast are initially assumed, 2) gravity are calculated for the assumed model and error between the true (observed) and calculated gravity are determined, 3) new vertices and density contrast are determined from the error by using the damped least-squares inversion method, and 4) final model is determined when the error is very small. Results of this study show that the shape and density contrast of each model are accurately determined when the density contrast is constant or vertical density gradient is known. In case where the density gradient is unknown, the inversion gives incorrect results. But the shape and density gradient of the model are determined when the surface density contrast is known.

• Journal of the Korean earth science society
• /
• v.23 no.1
• /
• pp.15-29
• /
• 2002

Bayesian inversion for gravity and resistivity data was performed to investigate the cavity structure appearing as a lava tunnel in Cheju Island, Korea. Dipole-dipole DC resistivity data were proposed for a prior information of gravity data and we applied the geostatistical techniques such as kriging and simulation algorithms to provide a prior model information and covariance matrix in data domain. The inverted resistivity section gave the indicator variogram modeling for each threshold and it provided spatial uncertainty to give a prior PDF by sequential indicator simulations. We also presented a more objective way to make data covariance matrix that reflects the state of the achieved field data by geostatistical technique, cross-validation. Then Gaussian approximation was adopted for the inference of characteristics of the marginal distributions of model parameters and Broyden update for simple calculation of sensitivity matrix and SVD was applied. Generally cavity investigation by geophysical exploration is difficult and success is hard to be achieved. However, this exotic multiple interpretations showed remarkable improvement and stability for interpretation when compared to data-fit alone results, and suggested the possibility of diverse application for Bayesian inversion in geophysical inverse problem.

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

Uiseong subbasin belonging to Kyungsang basin resulted from volcanic activity in the late Cretaceous. In this study, we carry out MT and gravity survey at the Hwasan caldera, which was formed of volcanic and abyssal rocks complex, then analyze and identify geological substructure. Potential survey such as gravity and magnetic survey has been mainly carried out in former studies, so depth information for understanding substructure was not enough. To complement a potential survey, we use MT method, which has high vertical resolution. Moreover we make a simple 2D model comparing with former study. The result of MT and gravity 2D modeling shows that this area is roughly composed of 3 layers; The bottom layer is a basement. In the second layer, intrusive rocks having high resistivity is placed along the ring faults and the sedimentary layer of low resistivity is inside caldera. The highest layer is alluvium. To comprehend the 3D structure of the Hwasan caldera, we perform 3D gravity inversion, and construct the 3D model from the result of 3D gravity inversion. MT responses are calculated by using the constructed 3D model and the 3D model of the Hwasan caldera's structure is suggested after comparing the calculated values with the observed values at MT line.

• Journal of the Korean earth science society
• /
• v.30 no.1
• /
• pp.49-57
• /
• 2009

We studied a structure of the Kunsan basin in the Yellow Sea using ship-borne magnetic data and altimetry satellite-derived gravity data provided from the Scripps institution of oceanography in 2006. The gravity data was analyzed via power spectrum analysis and gravity inversion, and the magnetic data via analytic signal technique, pseudo-gravity transformation, and its inversion. The results showed that the depth of bedrock tended to increase as we approached the center of the South Central Sag in Kunsan basin and that the maximum and minimum of its depth were estimated to be about 6-8 km and 2 km, respectively. Inaddition, the observed high anomaly of gravity and magnetism was attributed to the intrusion of igneous rock of higher density than the surrounding basement rock in the center of South Central Sag, which was consistent with the interpretation of seismic data obtained in the same region.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.7
• /
• pp.228-235
• /
• 2019

The purpose of this study was to investigate the effect of 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise on VAS, ROM, isokinetic myofunction, and dynamic stability after surgery of modified brostrom operation in chronic ankle instability patients. The subjects of this study were 12 chronic ankle instability patients who underwent modified Brostrom operation(MBO) by the same doctor. 6 weeks' accelerated rehabilitation program is scheduled to perform for 60min, everyday, and also anti-gravity program performed for 15~30min, everyday. The visual analog scale(VSA) and significantly decreased(p<.001) and ROM in all of dorsal flexion, plantar flexion, inversion and eversion significantly increased(p<.05) after 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise. Both inversion and eversion peak torque at $60^{\circ}/sec$(p<.001, p<.01) and at $180^{\circ}/sec$(p<.001) significantly increased after 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise respectively. In muscle defect, although inversion(p<.01) and eversion(p<.001) at $60^{\circ}/sec$ and inversion(p<.01) at $180^{\circ}/sec$ significantly decreased, eversion at $180^{\circ}/sec$ tended to decrease but did not change significantly after 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise. The dynamic stability significantly increased after 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise(p<.001). These results suggest that 6 weeks' accelerated rehabilitation with anti-gravity treadmill exercise has positive effect of VAS, ROM, isokinetic myofunction, and dynamic stability after surgery of modified brostrom operation in chronic ankle instability patients. Therefore, we consider that the accelerated rehabilitation with anti-gravity treadmill exercise, which is safely and fast method, has effect on more faster recovery of ankle stability, play ground and normal daily activities.

• Economic and Environmental Geology
• /
• v.28 no.4
• /
• pp.395-404
• /
• 1995

The geologic structure of the Cheju volcanic island has been investigated by analyzing the gravity and magnetic data. Bouguer gravity map shows apparent circular low anomalies at the central volacanic edifice, and the maximum difference of the anomaly values on the island appears to be 30 mgal. The subsurface structure of the island is modeled by three-dimensional depth inversion of gravity data by assuming the model consists of a stacked grid of rectangular prisms of volcanic rocks bounded below by basement rocks. The gravity modeling reveals that the interface between upper volvanic rocks and underlying basement warps downward under Mt. Halla with the maximum depth of 5 km. Magnetic data involve aeromagnetic and surface magnetic survey data. Both magnetic anomaly maps show characteristic features which resemble the typical pattern of total magnetic anomalies caused by a magnetic body magnetized in the direction of the geomagnetic field in the middle latitude region, though details of two maps are somewhat different. The reduced-to-pole magnetic anomaly maps reveal that main magnetic sources in the island are rift zones and the Halla volcanic edifice. The apparent magnetic boundaries inferred by the method of Cordell and Grauch (1985) are relatively well matched with known geologic boundaries such as that of Pyosunri basalt and Sihungri basalt which form the latest erupted masses. Inversion of aeromagnetic data was conducted with two variables: depth and susceptibility. The inversion results show high susceptibility bodies in rift zones along the long axis of the island, and at the central volcano. Depths to the basement are 1.5~3 km under the major axis, 1~1.5 km under the lava plateau and culminates at about 5 km under Mt. Halla. The prominent anomalies showing N-S trending appear in the eastern part of both gravity and magnetic maps. It is speculated that this trend may be associated with an undefined fault developed across the rift zones.