• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.171-182
• /
• 2018

We performed the seismic field data processing using an open-source software (Madagascar) to verify if it is applicable to processing of field data, which has low signal-to-noise ratio and high uncertainties in velocities. The Madagascar, based on Python, is usually supposed to be better in the development of processing technologies due to its capabilities of multidimensional data analysis and reproducibility. However, this open-source software has not been widely used so far for field data processing because of complicated interfaces and data structure system. To verify the effectiveness of the Madagascar software on field data, we applied it to a typical seismic data processing flow including data loading, geometry build-up, F-K filter, predictive deconvolution, velocity analysis, normal moveout correction, stack, and migration. The field data for the test were acquired in Gunsan Basin, Yellow Sea using a streamer consisting of 480 channels and 4 arrays of air-guns. The results at all processing step are compared with those processed with Landmark's ProMAX (SeisSpace R5000) which is a commercial processing software. Madagascar shows relatively high efficiencies in data IO and management as well as reproducibility. Additionally, it shows quick and exact calculations in some automated procedures such as stacking velocity analysis. There were no remarkable differences in the results after applying the signal enhancement flows of both software. For the deeper part of the substructure image, however, the commercial software shows better results than the open-source software. This is simply because the commercial software has various flows for de-multiple and provides interactive processing environments for delicate processing works compared to Madagascar. Considering that many researchers around the world are developing various data processing algorithms for Madagascar, we can expect that the open-source software such as Madagascar can be widely used for commercial-level processing with the strength of expandability, cost effectiveness and reproducibility.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.280-290
• /
• 2002

Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.223-235
• /
• 2002

Since weak Bones or geological lineaments are likely to be eroded, there may develop weak Bones beneath rivers, and a careful evaluation of ground condition is important to construct structures passing through a river. DC resistivity method, however, has seldomly applied to the investigation of water-covered area, possibly because of difficulties in data aquisition and interpretation. The data aquisition having high quality may be the most important factor, and is more difficult than that in land survey, due to the water layer overlying the underground structure to be imaged. Through the numerical modeling and the analysis of a case history, we studied the method of resistivity survey at the water-covered area, starting from the characteristics of measured data, via data acquisition method, to the interpretation method. We unfolded our discussion according to the installed locations of electrodes, ie., floating them on the water surface, and installing them at the water bottom, because the methods of data acquisition and interpretation vary depending on the electrode location. Through this study, we could confirm that the DC resistivity method can provide fairly reasonable subsurface images. It was also shown that installing electrodes at the water bottom can give the subsurface image with much higher resolution than floating them on the water surface. Since the data acquired at the water-covered area have much lower sensitivity to the underground structure than those at the land, and can be contaminated by the higher noise, such as streaming potential, it would be very important to select the acquisition method and electrode array being able to provide the higher signal-to-noise ratio (S/N ratio) data as well as the high resolving power. Some of the modified electrode arrays can provide the data having reasonably high S/N ratio and need not to install remote electrode(s), and thus, they may be suitable to the resistivity survey at the water-covered area.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.39-46
• /
• 2012

Recently, stratigraphic reservoirs are getting more attention than structural reservoirs which have mostly developed. However, recognizing stratigraphic thin gas reservoirs in a stacked section is usually difficult because of tuning effects. Moreover, if the reflections from the brine-saturated region of a thin layer have the same polarity with those from the gas-saturated region, we could not easily identify the gas reservoir with conventional data processing technique. In this study, we introduced a way to delineate the gas-saturated region in a thin layer reservoir using a spectral decomposition method. First of all, amplitude spectrum with the variation of the frequency and the incident angle was investigated for the medium which represents property of Class 3, Class 1 or Class 4 AVO response. The results show that the maximum difference in the amplitude spectra between brine and gas-saturated thin layers occurs around the peak frequency independent of the incident angle and the type of AVO responses. In addition, the amplitude spectra of the gas-saturated zone are greater than those of brine-saturated one in Class 3 and Class 4 at the peak frequency while those of phenomenon occur oppositely in Class 1. Based on the results, we applied spectral decomposition method to the stacked section in order to distinguish the gas-saturated zone from the brine-saturated zone in a thin layer reservoir. To verify our new method, we constructed a thin-layer velocity model which contains both gas and brine-saturated zones which have the same reflection polarities. As a result, in the spectral decomposed sections near the peak frequency obtained by Wigner-Ville Distribution (WVD), we could identify the difference between reflections from gas- and brinesaturated region in the thin layer reservoir, which was hardly distinguishable in the stacked section.

• The Korean Journal of Petroleum Geology
• /
• v.2 no.2 s.3
• /
• pp.59-70
• /
• 1994

The Cenozoic geological structures and the tectonic evolution of the southern Ulleung Basin were studied with seismic profiles and exploration well data. Basement structure of the Korea Strait is distinctly characterized by normal faults trending northeast to southwest. The normal faults of the basement are most likely related to the initial liking and extensional tectonics of Ulleung Basin. Tsushima fault along the west coast of Tsushima islands runs northeastward to the central Ulleung Basin. The Middle Miocene and older sequences in the Tsushima Strait show folds and faults mostly trending northeast to southwest. These folds and faults may be interpreted as a result of compressional tectonics. The Late Miocene to Qauternary sequences are not much deformed, but numerous faults mostly N-S trending are dominated in the Tsushima Strait. The Ulleung Basin was in intial rifting during Oligocene, and then active extension and subsidence from Early to early Middle Miocene. Therefore SW Japan separated from Korea Peninsula and drifted toward southeast, and Ulleung Basin was formed as a pull-apart basin under dextral transtensional tectonic regime. During rifting and extensional stage, Tsushima fault as a main tectonic line separating SW Japan block from the Korean Peninsula acted as a normal faulting with right-lateral strike-slip motion as SW Japan drifted southeastward. During middle Middle Miocene to early Late Miocene, the opening of Ulleung basin stopped and uplifted due to compressional tectonics. The southwest Japan block converging on the Korean Peninsula caused compressional stress to the southern margin of Ulleung Basin, resulting in strong deformation under sinistral transpressional tectonic regime. Tsushima fault acted as thrust fault with left-lateral strike-slip motion. From middle Late Miocene to Quaternary, the southern margin of Ulleung Basin has been controlled by compressional motion. Thus the Tsushima fault still appears to be an active thrust fault by compressional tectonic regime.

• The Korean Journal of Petroleum Geology
• /
• v.2 no.2 s.3
• /
• pp.51-57
• /
• 1994

Analysis of the aeromagnetic data aquired by US Navy in the year 1969 permits us to predict a new sedimentary basin, Heugsan Basin, south of the known Gunsan Basin in Block Ⅱ. The basin appears to consist of three sub-basins trending NNW-SSE. The results of our analysis provide not only an independent assessment of the Gunsan Basin, but also new important information on the tectonic origin and mechanism for the two basins as well as for the entire region. The basin forming tectonic style is interpreted as rhombochasm associated with double overstepped left-lateral wrench faults. From the magnetic evidence, a few NE-SW trending major onshore faults are extended to the study area. We also interpreted the nature of the faults to be left-lateral wrenches. This new gross structural style is consistent with the results of recent Yeongdong Basin analysis by Lee. The senses of fault movement are also supported by the paleomagnetic evidence that the Philippine Sea had experienced an 80-degree clockwise rotation since the Eocene. Based on a 2 $\frac{1}{2}$ model study the probable maximum thickness of the sediments in the Gunsan Basin is approximately 7500 meters. We believe that the new Heugsan Basin was left unidentified because a high velocity layer may be overlying the basin. Because the overall structural configuration of the Heugsan Basin appears to be favorable for hydrocarbon accumulation, a detailed airborne magnetic survey is recommended in the area in order to verify the magnetic expression of this thick basin. A detailed subsequent marine gravity survey is also recommended in order to delineate the sedimentary section and to acquire supplemental data to the magnetic method only if an overlying high velocity layer is confirmed. Otherwise a high energy source seismic survey may be more effective.

• Korean Journal of Remote Sensing
• /
• v.25 no.1
• /
• pp.31-44
• /
• 2009

The objective of this study was to measure backscattering coefficients of paddy rice using a L-, C-, and X-band scatterometer system with full polarization and various angles during the rice growth period and to relate backscattering coefficients to rice growth parameters. Radar backscattering measurements of paddy rice field using multifrequency (L, C, and X) and full polarization were conducted at an experimental field located in National Academy of Agricultural Science (NAAS), Suwon, Korea. The scatterometer system consists of dual-polarimetric square horn antennas, HP8720D vector network analyzer ($20\;MHz{\sim}20\;GHz$), RF cables, and a personal computer that controls frequency, polarization and data storage. The backscattering coefficients were calculated by applying radar equation for the measured at incidence angles between $20^{\circ}$ and $60^{\circ}$ with $5^{\circ}$ interval for four polarization (HH, VV, HV, VH), respectively. We measured the temporal variations of backscattering coefficients of the rice crop at L-, C-, X-band during a rice growth period. In three bands, VV-polarized backscattering coefficients were higher than hh-polarized backscattering coefficients during rooting stage (mid-June) and HH-polarized backscattering coefficients were higher than VV-, HV/VH-polarized backscattering coefficients after panicle initiation stage (mid-July). Cross polarized backscattering coefficients in X-band increased towards the heading stage (mid-Aug) and thereafter saturated, again increased near the harvesting season. Backscattering coefficients of range at X-band were lower than that of L-, C-band. HH-, VV-polarized ${\sigma}^{\circ}$ steadily increased toward panicle initiation stage and thereafter decreased, and again increased near the harvesting season. We plotted the relationship between backscattering coefficients with L-, C-, X-band and rice growth parameters. Biomass was correlated with L-band hh-polarization at a large incident angle. LAI (Leaf Area Index) was highly correlated with C-band HH- and cross-polarizations. Grain weight was correlated with backscattering coefficients of X-band VV-polarization at a large incidence angle. X-band was sensitive to grain maturity during the post heading stage.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.25-30
• /
• 1998

Time-lapse crosswell seismic data, recorded before and after the cavity filling, showed that the filling increased the velocity at a known cavity zone in an old mine site in Inchon area. The seismic response depicted on the tomogram and in conjunction with the geologic data from drillings imply that the size of the cavity may be either small or filled by debris. In this study, I attempted to evaluate the filling effect by analyzing velocity measured from the time-lapse tomograms. The data acquired by a downhole airgun and 24-channel hydrophone system revealed that there exists measurable amounts of source statics. I presented a methodology to estimate the source statics. The procedure for this method is: 1) examine the source firing-time for each source, and remove the effect of irregular firing time, and 2) estimate the residual statics caused by inaccurate source positioning. This proposed multi-step inversion may reduce high frequency numerical noise and enhance the resolution at the zone of interest. The multi-step inversion with different starting models successfully shows the subtle velocity changes at the small cavity zone. The inversion procedure is: 1) conduct an inversion using regular sized cells, and generate an image of gross velocity structure by applying a 2-D median filter on the resulting tomogram, and 2) construct the starting velocity model by modifying the final velocity model from the first phase. The model was modified so that the zone of interest consists of small-sized grids. The final velocity model developed from the baseline survey was as a starting velocity model on the monitor inversion. Since we expected a velocity change only in the cavity zone, in the monitor inversion, we can significantly reduce the number of model parameters by fixing the model out-side the cavity zone equal to the baseline model.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.22-37
• /
• 2020

We analyze the comprehensive three-dimensional (3D) magnetic structure characteristics from the seafloor to the deep layer of the Tofua Arc (TA) 12 seamount in the Tonga Arc, Southwestern Pacific, using bathymetric and geomagnetic data, and magnetization vector inversion (MVI) results. The seafloor features surrounding TA 12 highlight a NW-SE-oriented elliptical caldera at the summit of the seamount, two small cones in the depressed area. A large-scale sea valley is present on the western flank of the seamount, extending from these caldera cones to the southwestern base of the seamount. TA 12 seamount exhibits a low magnetic anomaly in the caldera depression, whereas a high magnetic anomaly is observed surrounding the low magnetic anomaly across the caldera summit. It is therefore presumed that there may be a strong magnetic material distribution or magma intrusion in the caldera. The 3D MVI results show that the high anomaly zones are mainly present in the surrounding slopes of the seamount from the seafloor to the -3,000 m (below the seafloor) level, with the magnetic susceptibility intensity increasing as the seafloor level increases at the caldera depression. However, small high anomaly zones are present across the study area near the seafloor level. Therefore, we expect that the magma ascent in TA 12 seamount migrated from the flanks to the depression area. Furthermore, we assume that the complex magnetic distribution near the seafloor is due to the remnant magnetization.