• Journal of the Korean Geophysical Society
• /
• v.2 no.2
• /
• pp.91-99
• /
• 1999

A high-resolution seismic profile acquired across the northeastern boundary of the Pungam Basin, one of the Cretaceous sedimentary basins in Korea, has been interpreted to delineate subsurface geological structures across the basin boundary. We identified boundary faults and unconformity surfaces of the basin and divided sediment body into three seismic depositional units (Units I, II, and III from youngest to oldest). Inferred from fault geometry and type, northeastern part of the Pungam Basin has been formed by a strike-slip fault whereas the normal faults near the boundary were formed by transtensional movement along a fault zone. A 350-400 m thick sediment layer is overlying the Precambrian gneiss. Bedding planes of Unit III are dipping westward and are closely related to an anticline in the acoustic basement. Unit II is also tilted westward, suggesting that the eastern part of the fault zone was uplifted after deposition of lower part of the sedimentary body. Afterward, the uplifted sediment layers were eroded and transported to the western part of the basin. Chaotic reflection pattern of sedimentary Units II and III may suggest that strike-slip movement along the fault zone deformed basin-filled sediments.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.23-32
• /
• 2002

Capacitively-coupled resistivity (CCR) system is known to be very useful where galvanic contact to earth is impossible, such as the area covered with thick ice, snow, concrete or asphalt. This system injects current non-galvanically, i.e., capacitively to earth through line antenna and measures potential difference in a same manner. We derived geometric factor for two types of antenna configuration and presented the method of processing and converting the data obtained with CCR system suitable to conventional resistivity inversion analysis. The CCR system, however, has limitations on use at conductive area or electrically noisy area since it is very difficult to inject sufficient current to earth with this system as with conventional resistivity system. This causes low SM ratio when acquiring data with CCR system and great care must be taken in acquiring data with this system. Additionally the uniform contact between line antennas and earth is also crucial factor to obtain good S/N ratio data. The CCR method, however, enables one to perform continuous profiling over a survey line by dragging entire system and thus will be useful in rapid investigation of conductivity distribution in shallow subsurface.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.39-47
• /
• 2000

A great deal of computing time and a large computer memory are needed to solve wave equation in a large complex subsurface layers using the finite difference method. The computing time and memory can be reduced by decreasing the number of grid points per minimum wave length. However, the decrease of grids may cause numerical dispersion and poor accuracy. In this study we performed the grid dispersion analysis for several rotated finite difference operators, which was commonly used to reduce grids per wavelength with accuracy in order to determine the solution for the acoustic wave equation in frequency domain. The rotated finite difference operators were to be extended to 81, 121 and 169 difference stars and studied whether the minimum grids could be reduced to 2 or not. To obtain accuracy (numerical errors less than $1\%$) the following was required: more than 13 grids for conventional 5 point difference stars, 9 grids for 9 difference stars, 3 grids for 25 difference stars, and 2.7 grids for 49 difference stars. After grid dispersion analysis for the new rotated finite difference operators, more than 2.5 grids for 81 difference stars, 2.3 grids for 121 difference stars and 2.1 grids for 169 difference stars were needed. However, in the 169 difference stars, there was no solution because of oscillation of the dispersion curves in the group velocity curves. This indicated that the grids couldn't be reduced to 2 in the frequency acoustic wave equation. According to grid dispersion analysis for the determination of grid points, the more rotated finite difference operators, the fewer grid points. However, the more rotated finite difference operators that are used, the more complex the difference equation terms.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.48-52
• /
• 2000

We have calculated and analyzed the electromagnetic responses of buried conductive pipes due to a horizontal magnetic dipole source on the pound using a three-dimensional (3-D) finite element method to provide useful guidelines for designing electromagnetic pipe locator and for field operation of the system. For single buried pipe, the horizontal component and the horizontal difference of the vertical component of magnetic field show peaks above the pipe. When comparing the width of response curves of both cases around the peak, horizontal difference of vertical component of magnetic field shows much narrower peak, 2 times narrower at a half of maximum amplitude, than that of horizontal component of magnetic field. Accordingly, we can pinpoint the horizontal location of pipe on the ground more accurately by measuring the horizontal difference of vertical component of magnetic fold. Moreover, it will have a merit in determining the depth of pipe, because the equation for depth estimation is defined just above the pipe. When there are two buried pipes separated by two meters with each other, the response of horizontal difference of vertical component of magnetic field has two separate peaks each of which is located above the pipe whereas horizontal magnetic field response has only one peak above the pipe just below the transmitter. Thus, when there exist more than a buried pipe, measuring the horizontal difference of vertical magnetic field can effectively detect not only the pipe under transmitter but also adjacent ones. The width of response curves also indicates higher resolving ability of horizontal difference of vertical component of magnetic field.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.304-315
• /
• 2006

We have examined the applicability of f-k analysis to the GPR direct wave measurement for water content to characterize vadose zone condition. When the vadose zone consists of a dry surface layer over wet substratum, we obtained f-k spectra where most of the energy is bounded by the air and dry soil velocities. In this case, dry soil velocity was successfully estimated by using high frequency data. On the other hands, when wet soil overlies dry substratum, the f-k spectra show a contrasting response where most of the energy travels with the velocity bounded by dry and wet soil velocities. In this case, the radar waves are trapped and guided within wet soil layer, exhibiting velocity dispersion. By adopting modal propagation theory, we could formulae a simple inversion code to find two layer's dielectric constants as well as layer thickness. By inverting the velocity dispersion curve obtained from f-k spectra of synthetic modeling data, we could obtain good estimates of dielectric constants of each layer as well as first layer thickness. Moreover, we could obtain more accurate results by including the higher mode data. We expect this method will be useful to get the quantitative property of real subsurface when the field condition is similar.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.187-201
• /
• 2014

Compared with the separated inversion of electromagnetic (EM) and seismic data, a joint inversion using both EM and seismic data reduces the uncertainty and gives the opportunity to use the advantage of each data. Seismic fullwaveform inversion allows velocity information with high resolution in complicated subsurface. However, it is an indirect survey which finds the structure containing oil and gas. On the other hand, marine controlled-source EM (mCSEM) inversion can directly indicate the oil and gas using different EM properties of hydrocarbon with marine sediments and cap rocks whereas it has poor resolution than seismic method. In this paper, we have developed a joint EM inversion algorithm using a cross-gradient technique. P-wave velocity structure obtained by full-waveform inversion using plane wave encoding is used as structure constraints to calculate the cross-gradient term in the joint inversion. When the jointinversion algorithm is applied to the synthetic data which are simulated for subsea reservoir exploration, images have been significantly improved over those obtained from separate EM inversion. The results indicate that the developed joint inversion scheme can be applied for detecting reservoir and calculating the accurate oil and gas reserves.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.76-85
• /
• 1993

In order to predict the free surface signature of turbulent ship wakes two things are essential; a basic understanding of the mechanism of turbulent vortical flow/free surface interactions and a mathematical model to accurately predict the signature. The goal of the study described here is both to supplement experimental work to obtain basic understanding, as well as to condense this understanding in a model(or models) that captures the essential phenomena and thus allows predictions. To do so we followed two main paths guided by experimental observations. One is full simulations of the flow using the clavier-Stokes equations. The other is a vortex modeling, where the vortical structures of the flows are approximated by idealized structures, an the interaction assumed to be essentially inviscid. These approaches complement each other. Full simulations are only applicable to small scale phenomena, where the system is simple, and the Reynolds number is low. The vortex modeling, on the other hand, cannot represent essentially viscous aspects of the problem such as the effect of contamination gradient. Obviously, the modeling is what may eventually lead to a prediction method; the full simulations-too limited to mimic all but the simplest circumstances-are to aid and support the construction of realistic models. We address two-dimensional aspects of the vortex/free surface interaction first. Secondly we obtain some basic understanding of the interaction process through an experiment and then talk about several three-dimensional problems hoping to develop a successful prediction model.