• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.191-196
• /
• 2004

Seismic waveform inversion can be solved by using the classical Gauss-Newton method, which needs to construct the huge Hessian by the directly computed Jacobian. The property of Hessian mainly depends upon a source and receiver aperture, a velocity model, an illumination Bone and a frequency content of source wavelet. In this paper, we try to invert the Marmousi seismic data by controlling the huge Hessian appearing in the Gauss-Newton method. Wemake the two kinds of he approximate Hessian. One is the banded Hessian and the other is the approximate Hessian with automatic gain function. One is that the 1st updated velocity model from the banded Hessian is nearly the same of the result from the full approximate Hessian. The other is that the stability using the automatic gain function is more improved than that without automatic gain control.

• Economic and Environmental Geology
• /
• v.49 no.3
• /
• pp.225-242
• /
• 2016

Rock physics serves as a useful tool for seismic reservoir characterization and monitoring by providing quantitative relationships between rock properties and seismic data. Rock physics models can predict effective moduli for reservoirs with different mineral components and pore fluids from well-log data. The distribution of reservoirs and fluids for the entire seismic volume can also be estimated from rock physics models. The first part of this report discusses the Voigt, Reuss, and Hashin-Shtrikman bounds for effective elastic moduli and the Gassmann fluid substitution. The second part reviews various contact models for moderate- to high-porosity sands. In the third part, constant-cement model, known to work well for the sand that gradually loses porosity with deteriorating sorting, was applied to the well-log data from an oil field in the North Sea. Lastly, the rock physics template constructed from the constant-cement model and the results from the prestack inversion of 2D seismic data were combined to predict the lithology and fluid types for the sand reservoir of this oil field.

• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.49-57
• /
• 2014

3D magnetic inversion, based on the assemblage of 2D forward modeling and inversion as a practical technique to reflect the a priori information, was conducted to investigate the spatial distribution features of black-shale related and pyrometamorphic uranium deposit and several lithological units of Ogcheon Super Group in an area of Geumsan. By using the 3D visualization technique with suitable susceptibility interval and horizontal slice map, the spatial distribution of magnetic susceptibility corresponded to the black shale related uranium bearing lithological units, Black Slate member was well coincided with a information of uranium deposit. Also, even though it is indirect indicator for the detetction of uranium deposits interbedded in Gray Hornfels member, spatial susceptibility distribution which shows the south-east magnetic linearment corresponding to the Majeon-ri formation and Dark Gray Slate were matched well. From this investigation, we inferred that maximum depth extension which Black Slate member can be separately recognized with respect to adjacent Dark Gray Slate with strong magnetic susceptibility anomaly is about 150m with reference elevation level of 306m. In addition, Majeon-ri formation located south of Black Slate member revels relatively high magnetic susceptibility range but shows high spatial susceptibility fluctation. And, as an intrusive rocks, Jurassic Biotite Granite shows relatively low magnetic susceptibility characteristics. On the contrary, Cretaceous granite distributed in soutthern part of the study area shows the relatively high susceptibility distribution.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.286-301
• /
• 2008

When magnetotelluric (MT) data are obtained in vicinity of the coast, the sea can distort observed MT responses, especially those of deep part of subsurface. We introduce an iterative method to correct the sea effect, based on the previous topographic correction method which removes the distortions due to topographic changes in seafloor MT data. The method first corrects the sea effect in observed MT impedance, and then inverts corrected responses in a model space without the sea. Due to mutual coupling between sea and subsurface structure, the correction and inversion steps are iterated until changes in each result become negligible. The method is validated for 1-D and 2-D structure using synthetic MT data produced by 3-D forward modeling including surrounding seas. In all cases, the method closely recovers the given structure after a few iterations. To test the applicability of the proposed method to field data, we generate synthetic MT data for the Jeju Island whose 1-D conductivity structure is well known, using 3-D forward modeling. The distortions due to the surrounding sea start to appear below the frequency about 1 Hz, and are relatively severe in the electrical field perpendicular to the coastline because of the location of the observation sites. The proposed method successfully eliminates the sea effect after three iterations, and both 1-D and 2-D inversion of corrected responses closely recover the given subsurface structure of the Jeju Island model.

• Journal of the Korean Geophysical Society
• /
• v.3 no.4
• /
• pp.235-244
• /
• 2000

We studied basic ideas of a network based Gravity/Magnetic data processing server/client system which provides functions of data processing, forward modeling, inversion and data process on Data Base. This Java technology was used to provide facilities, socket communication and JDBC(Java Database Connectivity) technology to produce an effective and practical client application. The server computers are linked by network to process the MPI parallelized computing. This can provide useful devices of the geophysical process and modeling that usually require massive computing performance and time. Since this system can be accessed by lots of users, it can provides the consistent and confident results through the verified processing programs. This system also makes it possible to get results and outputs through internet when their local machines are connected to the network. It can help many users who want to omit the jobs of system administration and to process data during their field works.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Proceedings of the Korean Society of Surveying, Geodesy, Photogrammetry, and Cartography Conference
• /
• 2002.10a
• /
• pp.195-199
• /
• 2002

SAR(Synthetic Aperture Radar) 영상은 경사촬영을 수행하므로 지형의 기복에 따른 영향을 많이 받는다. 따라서 SAR영상을 이용하여 여러 가지 정보들을 추출하여 이용하기 위해서는 전처리 과정으로서 지형의 기복에 따른 여러 가지 왜곡들을 보정해야 한다. 이에 본 연구에서는 RADARSAT SAR 영상을 이용하여 궤도모델링, 정사보정을 수행하고 역산란계수, 국부입사각 계산 등을 통해 지형기복에 따른 방사왜곡보정을 수행하였다.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.250-262
• /
• 2008

In order to assess the feasibility of the dipole-dipole electric method to the investigation of metallic ore deposit, both field data simulation and inversion are carried out for several simplified ore deposit models. Our interest is in a vein-type model, because most of the ore deposits (more than 70%) exist in a vein type in Korea. Based on the fact that the width of the vein-type ore deposits ranges from tens of centimeters to 2m, we change the width and the material property of the vein, and we use 40m-electrode spacing for our test. For the vein-type model with too small width, the low resistivity zone is not detected, even though the resistivity of the vein amounts to 1/300 of that of the surrounding rock. Considering a wide electrode interval and cell size used in the inversion, it is natural that the size of the low resistivity zone is overestimated. We also perform field data simulation and inversion for a vein-type model with surrounding hydrothermal alteration zones, which is a typical structure in an epithermal ore deposits. In the model, the material properties are assumed on the basis of resistivity values directly observed in a mine originated from an epithermal ore deposits. From this simulation, we can also note that the high resistivity value of the vein does not affect the results when the width of the vein is narrow. This indicates that our main target should be surrounding hydrothermal alteration zones rather than veins in field survey. From these results, we can summarize that when the vein is placed at the deep part and the difference of resistivity values between the vein and the surrounding rock is not large enough, we cannot detect low resistivity zone and interpret the subsurface structures incorrectly using the electric method performed at the surface. Although this work is a little simple, it can be used as references for field survey design and field data Interpretation. If we perform field data simulation and inversion for a number of models and provide some references, they will be helpful in real field survey and interpretation.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.221-229
• /
• 2008

The three dimensional model method is widely applied in resource development for feasibility study, mine design, excavation planning and process management by constructing the database of various data in 3 dimensional space. Most of geophysical surveys for the purpose of engineering and resource development are performed in 2 dimensional line survey due to the restriction of the field situation, technical or economical situation and so on. The acquired geophysical data are used as the input for the 2 dimensional inversion under the 2 dimensional assumption. But the geophysical data are affected by 3 dimensional space. Therefore in order to reduce the error caused by 2 dimensional assumption, the 2 dimensional inversion result must be interpreted considering the additional information such as 3 dimensional topography, geological structure, borehole survey etc. The applicability and usability of 3 dimensional modeling method are studied by reviewing the case study to the geophysical data acquired in field of engineering and resource development.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2010.10a
• /
• pp.33-33
• /
• 2010

남극 반도 북부에 설치된 임시지진관측소와 세종기지 지진관측소에서 관측된 원격지진을 이용하여 P파의 상부 맨틀 속도 구조 모델링을 하였다. 사용된 자료는 1997년부터 1999년까지 실시된 SEPA (the Seismic Experiment in Patagonia and Antarctica) 연구에 설치된 7대의 임시 광대역 관측소와 IRIS/GSN 영구 관측소인 PMSA, 그리고 칠레 Jubany 기지(JUBA)와 아르헨티나의 Esperanza 기지(ESPZ)에 설치된 광대역 지진관측소에서 관측된 자료를 사용하였다. 모든 관측소가 남극 반도와 남 셔틀랜드 군도에 위치하고 있기 때문에 매우 낮은 신호대 잡음비를 보여주고 있다. 모델링에 사용된 자료는 95개 지진에서 축출한 347개의 P파와 PKP파로 실시되었다. 역산된 상부 맨틀의 속도 구조는 남 셔틀랜드 군도의 북쪽에서 빠르고 브랜스필드 스트레이 지역에서는 느린 속도 구조를 보여주고 있다.