• 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파로 실시되었다. 역산된 상부 맨틀의 속도 구조는 남 셔틀랜드 군도의 북쪽에서 빠르고 브랜스필드 스트레이 지역에서는 느린 속도 구조를 보여주고 있다.

• Journal of the Korean Geotechnical Society
• /
• v.23 no.2
• /
• pp.5-17
• /
• 2007

Various field setup and filtering criteria have been suggested to avoid the near field effects in surface wave methods. Unlike other surface wave methods HWAW method uses the near field component positively. It is possible by using maximum energy point based on time-frequency map and inversion method to consider receiver locations from the source point and body wave component. To verify the HWAW method in the near field numerical study was performed and the wave propagation in the stratified soil media was simulated due to a surface point load. All of five representative soil models were used. The experimental dispersion curves, determined by HWAW method at the various receiver distances in the region of near field, all coincided well with the theoretical dispersion curves determined by 3D forward modeling (Kausel's method). Consequently, it was considered that the HWAW method can provide reliable $V_s$ profiles effectively in the near field.

• Geophysics and Geophysical Exploration
• /
• v.14 no.2
• /
• pp.152-157
• /
• 2011

The small-loop electromagnetic (EM) method is one of the rapid and non-destructive geophysical methods and has been used widely for many geophysical investigations, particularly for shallow engineering and environmental surveys. Especially in the shallow marine environment, the small-loop EM technique is very effective because of rapid and convenient data acquisition, large signal and low noise level. However, the method has been rarely applied in the very conductive marine environment since it's penetration or investigation depth might be considered too low. In this study, we demonstrated that the small-loop EM method can be effectively applied in the extremely conductive marine environment through the analysis of 1D small-loop EM data. Furthermore, we confirmed that the resistivity distribution under the sea bottom can be quantitatively predicted from the 1D inversion results of synthetic and field data.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.327-335
• /
• 2006

Three processing strategies of seismic refraction data are tested in terms of velocity and depth profiles or structures for mapping of geological discontinuities: GRM(generalized reciprocal method), GLI(generalized linear inversion), Tomography. The test data used in this study are the shot gathers reconstructed by numerical modeling for the structures of 3 planar layers(horizontal, inclined), the buried vertical fracture zones, and vertical fault zones. Tomography is shown to be very efficient for mapping of more complicated tone such as vertical fault and buried fracture zones, whereas GRM and GLI can be useful for horizontal and/or inclined layers, probably on their bases of analysis of first arrivals in travel time curves.

• The Journal of Engineering Geology
• /
• v.17 no.2 s.52
• /
• pp.271-277
• /
• 2007

To find out the anomalous zone in cut slope composed of phyllite and shist, we performed resistivity tomography using a pole-dipole way. The electrical distribution that propagates from a current source in lower part of slope is measured by a potential electrode in upper part of slope. Apparent resistivity data are inverted with an iterative regularized inversion method to reconstruct 3D resistivity image. By comparing with the resistivity images in relation to each section, the images of anomalous zone correspond to their positions represented in cut slope. Therefore, the application of resistivity tomography in cut slope is useful to recognize the extension of anomalous zone.

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

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

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2009.10a
• /
• pp.75-80
• /
• 2009

Despite its flexibility to complex geometry, three-dimensional (3D) electromagnetic(EM) modeling schemes using finite element method (FEM) have been faced to practical limitation due to the resulting large system of equations to be solved. An efficient 3D FEM modeling scheme has been developed, which can adopt either direct or iterative solver depending on the problems. The direct solver PARDISO can reduce the computing time remarkably by incorporating parallel computing on multi-core processor systems, which is appropriate for single frequency multi-source configurations. When limited memory, the iterative solver BiCGSTAB(1) can provide fast and stable convergence. Efficient 3D simulations can be performed by choosing an optimum solver depending on the computing environment and the problems to be solved. This modeling includes various types of controlled-sources and can be exploited as an efficient engine for 3D inversion.