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

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.102-115
• /
• 2012

The purpose of this paper is to review several different methods calculating gas hydrate saturations. There are three methods using downhole log data, core data (including pressure core), and seismic velocity data. Archie's equation using electrical resistivity of downhole log data is widely used for saturation calculation. In this case, Archie's parameters should be defined accurately. And the occurrence types of gas hydrate significantly affect to saturation calculation. Thus saturation calculation should be carefully conducted. The methods using chlorinity and pressure core data are directly calculated from core sample. So far, the saturation calculated from pressure core gives accurate and quantitative values. But this method is needed much more time and cost. Thus acquisition of the continuous data with sediment depth is realistically hard. The recent several results show that the saturation calculated from resistivity data is the highest values, while the value calculated from pressure core is the lowest. But this trend is not always absolutely. Thus, to estimate accurate gas hydrate saturation, the values calculated from several methods should be compared.

• Geophysics and Geophysical Exploration
• /
• v.23 no.4
• /
• pp.243-252
• /
• 2020

To share an understanding of trajectory measurement in surveys using borehole, this tutorial summarizes the relevant mathematical principles of the borehole deviation survey based on coordinate transform. For uncased or open holes, calculations of the azimuth-deviation-tool face rotation using three-component accelerometer and magnetometer measurements are summarized. For the steel-cased holes, calculations are based on the time-derivative formula of the coordinate transform matrix; yaw-pitch-roll angles through time are mathematically determined by integrating the threecomponent angular velocity measurements from the gyroscope while also removing the Earth's rotation effect. Sensor and data fusion to increase the accuracy of borehole deviation survey is explained with an example of the method. These principles of borehole deviation surveys can be adapted for attitude estimation in air-borne surveys or for positioning in tunnels where global positioning system (GPS) signals cannot be accessed. Information on the optimization filter that must be incorporated in sensor fusion is introduced to help future research.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.17-25
• /
• 1999

A traveltime tomography has been carried out by transforming electromagnetic data in frequency domain to wave-like domain. The transform uniquely relates a field satisfying a diffusion equation to an integral of the corresponding wavefield. But direct transform of frequency domain magnetic fields to wave-field domain is ill-posed problem because the kernel of the integral transform is highly damped. In this study, instead of solving such an unstable problem, it is assumed that wave-fields in transformed domain can be approximated by sum of ray series. And for further simplicity, reflection and refraction energy compared to that of direct wave is weak enough to be neglected. Then first arrival can be approximated by calculating the traveltime of direct wave only. But these assumptions are valid when the conductivity contrast between background medium and the target anomalous body is low enough. So this approach can only be applied to the models with low conductivity contrast. To verify the algorithm, traveltime calculated by this approach was compared to that of direct transform method and exact traveltime, calculated analytically, for homogeneous whole space. The error in first arrival picked by this study was less than that of direct transformation method, especially when the number of frequency samples is less than 10, or when the data are noisy. Layered earth model with varying conductivity contrasts and inclined dyke model have been successfully imaged by applying nonlinear traveltime tomography in 30 iterations within three CPU minutes on a IBM Pentium Pro 200 MHz.

• Journal of the Mineralogical Society of Korea
• /
• v.19 no.4 s.50
• /
• pp.221-229
• /
• 2006

Mineralogical and chemical examinations were performed on interstratified illite-smectite (I-S) minerals that occur in the mudstones from a petroleum exploration well in the Tertiary marine basin, Japan. X-ray diffraction analysis shows that component layers of illite in the interstratified I-S increase with increasing burial depth while those of smectie decrease. In addition, the randomly (R=0) interstratified illite-smectite is changed into Rp1 ordered I-S at a depth of about 4,000 m, which corresponds to the result of organic analysis and indicates a burial temperature of about $100^{\circ}C$. However, the present geothermal gradient shows that the conversion of the random I-S to R=0 ordered I-S is likely to occur at 3,000 m. This discrepancy may be interpreted by the reverse fault at 2,500 m which resulted in a deeper burial of sediments up to 1,000 m. Chemical analysis also shows the compositional variation in I-S with increasing depth: a decrease in Si and an increases in Al and K, indicating that the substitution of Al for Si in tetrahedral sheets is compensated by the addition of K to interlayers. K may be derived from K-feldspar and micas, which is present in the mudstones.

• Geophysics and Geophysical Exploration
• /
• v.18 no.1
• /
• pp.31-34
• /
• 2015

A gravity meter has been used for exploring subsurface mineral resources and monitoring long-period events such as Earth tides. Recently, researchers found several other intriguing features that we could even detect large teleseismic earthquakes and monitor seismic ambient noise using gravimeters. The zero-length spring suspension technology gives the gPhone (Micro-g LaCoste) excellent low frequency sensitivity, which may have implications for investigating much longer-period natural events (e.g., Earth's hum, tsunami waves, etc.). In this study, we present preliminary results through temporary operation of the gPhone at Geumsan in South Korea for 9 months (Nov. 2008-Jul. 2009). The gPhone successfully recorded large teleseismic events and showed a clear seasonal variation of the Double frequency microseisms during its operation period.

• Korean Journal of Remote Sensing
• /
• v.33 no.3
• /
• pp.313-324
• /
• 2017

A lava tube is one of the hot issues of lunar science because it is regarded as a good candidate place for setting a lunar base. Recently much effort has been made to find lunar lava tubes. However, preceding works mainly made use of high-resolution lunar surface image data in conjunction with geomorphological consideration to present some lava tube candidates. Yet, those candidates stay no more than indirect indications. We propose a new data analysis technique of High Frequency (HF) radar observation data to find lunar lava tubes of which location depth is smaller than the range resolution of the radar pulse. Such shallow target echoes cannot be resolved from surface echoes, which presents the different location of the lunar surface compared to that of real lunar surface. The proposed technique instead finds the surface range (distance from LRS to the reflector of the most intense signal) anomaly which occurs as a result of the low range resolution of LRS pulse. We applied this technique to the surface range of Kaguya Lunar Radar Sounder (LRS) data. The surface range was deduced to make LRS surface elevation which was compared with the average surface elevation of Kaguya Digital Terrain Model (DTM). An anomalous discrepancy of the surface elevation was found in the Rima Galilaei area, which suggests the existence of a shallow lava tube.

• Journal of Ocean Engineering and Technology
• /
• v.30 no.5
• /
• pp.389-399
• /
• 2016

This paper presents the explorations of hydrothermal vents located in the Marina Arc and Back Arc Basin using the deep-sea ROV Hemire. These explorations were conducted by KRISO and KIOST to demonstrate the capability of Hemire in various applications for deep-sea scientific research. The missions included the following: (1) to search the reported vents, (2) conduct visual inspections, (3) deploy/recover a sediment trap and bait traps, (4) sample sediment/water/rock, (5) measure the magnetic field at the vent site, and (6) acquire a detailed map using multi-beam sonar near the bottom. We installed three HD cameras for precise visual inspection, a high-temperature thermometer, a three-component magnetometer, and a multi-beam sonar to acquire details of the bottom contour or identify vents in the survey area. The explorations were performed in an expedition from March 23 to April 5, 2016, and the missions were successfully completed. This paper discusses the operational process, navigation, and control of Hemire, as well as the exploration results.

• Geophysics and Geophysical Exploration
• /
• v.22 no.3
• /
• pp.116-131
• /
• 2019

In order to obtain high-quality high-resolution marine seismic data, the survey needs to be carried out at very low-sea condition. However, the survey is often performed with a slight wave, which degrades the quality of data. In this case, it is possible to improve the quality of seismic data by detecting the exact location of the sea bottom signal and eliminating the influence of waves or swells automatically during data processing. However, if noise is included or the sea bottom signal is weakened due to sea waves, sea bottom detection errors are likely to occur. In this study, we applied a method reducing such errors by estimating the sea bottom location, setting a narrow detection range and detecting the sea bottom location within this range. The expected location of the sea bottom was calculated using previously detected sea bottom locations for each channel of multi-channel data. The expected location calculated in each channel is also compared and verified with expected locations of other channels in a shot gather. As a result of applying this method to the noisy 8-channel high-resolution air-gun seismic data acquired off Yeosu, the errors in selecting the strong noise before sea bottom or the strong subsurface reflected signal after the sea bottom signal are remarkably reduced and it is possible to produce the high-quality seismic section with the correction of ~ 2.5 m swell effect.

• Geophysics and Geophysical Exploration
• /
• v.25 no.4
• /
• pp.203-213
• /
• 2022

To understand an underlying geological structure via seismic imaging, the velocity information of the subsurface medium is crucial. Although the full-waveform inversion (FWI) method is considered useful for estimating subsurface velocity models, 3D FWI needs a lot-of computing power and time. Herein, we compare the calculation efficiency and accuracy of frequency-domain 2D and 3D acoustic FWIs. Thereafter, we demonstrate that the artifacts from 2D approximation can be partially suppressed via frequency-domain 2D FWI by employing diffraction angle filtering (DAF). By applying DAF, which employs only big reflection angle components, the impact of noise and out-of-plane reflections can be reduced. Additionally, it is anticipated that the DAF can create long-wavelength velocity structures for 3D FWI and migration.