• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.315-322
• /
• 2010

The phase-screen and the split-step Fourier migrations, which are implemented in both the frequency-wavenumber and frequency-space domains by using one-way scalar wave equation, allow imaging in laterally heterogeneous media with less computing time and efficiency. The generalized-screen migration employs the series expansion of the exponential, unlike the phase-screen and the split-step Fourier migrations which assume the vertical propagation in frequency-wavenumber domain. In addition, since the generalized-screen migration generalizes the series expansion of the vertical slowness, it can utilize higher-order terms of that series expansion. As a result, the generalized-screen migration has higher accuracy in computing the propagation with wide angles than the phase-screen and split-step Fourier migrations for media with large and rapid lateral velocity variations. In this study, we developed a 2D prestack generalized-screen migration module for imaging a complex subsurface efficiently, which includes various dips and large lateral variations. We compared the generalized-screen propagator with the phase-screen propagator for a constant perturbation model and the SEG/EAGE salt dome model. The generalized-screen propagator was more accurate than the phase-screen propagator in computing the propagation with wide angles. Furthermore, the more the higher-order terms were added for the generalized-screen propagator, the more the accuracy was increased. Finally, we compared the results of the generalizedscreen migration with those of the phase-screen migration for a model which included various dips and large lateral velocity variations and the synthetic data of the SEG/EAGE salt dome model. In the generalized-screen migration section, reflectors were positioned more accurately than in the phase-screen migration section.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.307-314
• /
• 2010

Although conventional seismic data processing is based on the assumption that the media are isotropic, the subsurface is often anisotropy in shale formation or carbonate with cracks and fractures. This paper presents the anisotropic parameter and seismic modeling in transversely isotropic media with a vertical symmetry axis using seismic physical modeling. The experiment was successfully carried out with VTI media, laminated bakelite material, using contact transducer of p and s-wave transmission. The variation of velocities with angle of incidence was clearly shown in anisotropic material. Comparing these velocities with the calculated phase velocities, the (P) and (S)-wave velocity observed in anisotropic material was a very good agreement with the calculated values. Anisotropic parameter ${\varepsilon}$, ${\delta}$, ${\gamma}$ was estimated by using Lame's constant calculated from the observed velocity. For the purpose of testing (S)-wave polarization, a birefringence experiment was carried out. The higher velocity was associated with the polarization parallel to the fracture, and the lower velocity was associated with the polarization perpendicular to the fracture.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.332-344
• /
• 2007

Various geophysical well logs have been made along the four deep wells in Pohang, Gyeongbuk. The primary focus of geophysical well loggings was to improve understanding the subsurface geologic structure, to evaluate in situ physical properties, and to estimate aquifer production zones using fluid temperature and conductivity gradient logs. Especially natural gamma logs interpreted with core logs of borehole BH-1 were useful to discriminate the lithology and to determine the lithologic sequences and boundaries consisting of semi-consolidated Tertiary sediments and intrusive rocks such as basic dyke and Cretaceous sediments. Cross-plot of physical properties inferred from geophysical well logs were used to identify rock types such as Cretaceous sandstone and mudstone, Tertiary sediments, rhyolite, and basic dyke. The temperature log indicated $82.51^{\circ}C$ at the depth of 1,981.3 meters in borehole BH-4. However, considering the temperature of borehole BH-2 measured under stable condition, we expect the temperature at the depth in borehole BH-4, if it is measured in stable condition, to be about 5 or $6^{\circ}C$ higher. Several permeable fractures also have been identified from temperature and conductivity gradient logs, and cutting logs.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.109-115
• /
• 2008

It is very important to estimate the physical properties of survey area and delineate the geological basement in marine site survey for the design of offshore structures. For the purpose of providing high quality data by means of engineering site survey, it is necessary to apply several survey techniques and carry out the integrated interpretation to each other. In this study, we applied single channel seismic reflection method and OBC (Ocean Bottom Cable) type seismic refraction method at shallow marine. We used a dual boomer-single channel streamer as a source-receiver in seismic reflection survey and airgun source-the developed OBC type streamer in seismic refraction survey. We made 24 channels OBC type streamer which has 4m channel interval and each channel is composed of single hydrophone and preamplifier. We tested the field applicability of the proposed method and applied the typical seismic data processing methods to the obtained reflection data in order to enhance the data quality and image resolution. In order to estimate the geological velocity distribution from refraction data, seismic refraction tomography technique was applied. Therefore, we could successfully perform time-depth conversion using the velocity information as an integrated interpretation. The proposed method could provide reliable geologic information such as sediment layer thickness and 3D basement depth map.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.153-160
• /
• 2008

Seismic modeling is used to simulate wave propagation in the earth. Although the earth's subsurface is usually semi-infinite, we cannot handle the semi-infinite model in seismic modeling because of limited computational resources. For this reason, we usually assume a finite-sized model in seismic modeling. In that case, we need to eliminate the edge reflections arising from the artificial boundaries introducing a proper boundary condition. In this study, we changed three kinds of boundary conditions (sponge boundary condition, Clayton and Engquist's absorbing boundary condition, and Higdon's transparent boundary condition) so that they can be applied in elastic wave modeling for anisotropic media. We then apply them to several models whose Poisson's ratios are different. Clayton and Engquist's absorbing boundary condition is unstable in both isotropic and anisotropic media, when Poisson's ratio is large. This indicates that the absorbing boundary condition can be applied in anisotropic media restrictively. Although the sponge boundary condition yields good results for both isotropic and anisotropic media, it requires too much computational memory and time. On the other hand, Higdon's transparent boundary condition is not only inexpensive, but also reduce reflections over a wide range of incident angles. We think that Higdon's transparent boundary condition can be a method of choice for anisotropic media, where Poisson's ratio is large.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.268-274
• /
• 2007

Traditional two-dimensional (2D) interpretation of magnetotelluric (MT) data utilizes only transverse magnetic (TM)-mode data, because 2D inversion of transverse electric (TE)-mode data results in spurious features when 3D structures exist in the subsurface. The application of a 3D inversion algorithm to a single MT profile can reduce contamination due to off-profile anomalies and help us to incorporate TE-mode data in the interpretation. In this study, we conduct 2D and 3D inversions of MT data observed along two lines in Jeju Island. First, we invert apparent resistivities and phases in the TM and TE modes separately. Then, we perform 2D joint inversion of both TM- and TE-mode data and 3D inversion of both Zxy- and Zyx-mode data corresponding to TE- and TM-mode data in 2D. The resistivity images derived from all four data show that the geoelectrical structure in Jeju Island is a three-layered earth with the resistive-conductive-resistive stratigraphy within a depth of 5 km. The 3D inversion does not produce clear anomalies in the reconstructed profile image, while all of 2D do. This attributed to the possibility that 2D inversion results are distorted by exiting off-profile 3D anomalies in Jeju. With 3D inversion of 2D profile MT data, we can deduce more reliable results that are not seriously distorted by off-profile 3D anomalies.

• Geophysics and Geophysical Exploration
• /
• v.11 no.2
• /
• pp.85-92
• /
• 2008

Magnetotelluric (MT) methods are widely applied as an effective exploration technique to geothermal surveys. Two-dimensional (2-D) analysis is frequently used to investigate a complicated subsurface structure in a geothermal region. A 2-D finite-element method (FEM) is usually applied to the MT analysis, but we must pay attention to the accuracy of so-called auxiliary fields. Rodi (1976) proposed an algorithm of improving the accuracy of auxiliary fields, and named it as the MOM method. Because it introduces zeros into the diagonal elements of coefficient matrix of the FEM total equation, a pivoting procedure applied to the symmetrical band matrix makes the numerical solution far less efficient. The MOM method was devised mainly for the inversion analysis, in which partial derivatives of both electric and magnetic fields with respect to model parameters are required. In the case of forward modeling, however, we do not have to resort to the MOM method; there is no need of modifying the coefficient matrix, and the auxiliary fields can be elicited from the regular FEM solution. The computational efficiency of the MOM method, however, can be greatly improved through a sophisticated rearrangement of the total equation.

• Proceedings of the International Union of Geodesy And Geophysics Korea Journal of Geophysical Research Conference
• /
• 2003.05a
• /
• pp.20-20
• /
• 2003

In the summers of 1997 and 1998 and in February of 2000 we made 570 measurements of the ambient geomagnetic field 120 cm above the pavement surface of State Route 130, south of Pahoa, the island of Hawaii using a three-component fluxgate magnetometer. We measured at every 15.2 m (50 feet) interval covering a distance of 6, 310 m (20, 704 ft) where both historic and pre-historic highly magnetic basalt flows underlie. We also collected 197 core samples from eight road cuts, 489 specimens of which were subject to AF demagnetizations at 5 - 10 mT level up to a maximum field of 60 mT. We observed significant inclination anomalies ranging from a minimum of $31^{\circ}$ to a maximum $40^{\circ}$ where a uniform inclination value of $36.7^{\circ}$ (International Geomagnetic Reference Field, IGRF) was expected. Since the mean of the observed inclinations is approximately $35^{\circ}$ we assume that the study area is slightly affected by the magnetic terrain effect to a systematically shallower inclinations for being located in the regionally sloping surface of the southern side of the island (Baag, et al., 1995). We observed inclination anomalies showing wider (spacial) wavelength (160 - 600 m) and higher amplitudes in the historic lava flows area than in the northern pre-historic flows. Our observations imply that preexisting inclination anomalies such as those that we observed would have been interpreted as paleosecular variation (PSV). These inclination anomalies can best be attributed to concealed underground highly magnetic dikes, channel type lava flows, on-and-off hydrothermal activities through fissure-like openings, etc. Both the within- and between-site dispersions of natural remanent magnetization (NRM) are largest (up to ${\pm}7^{\circ}$) above the flows of 1955, while the area of pre-historic flows in the northern part of the study area exhibit the smallest dispersion. Nevertheless, mean inclinations of each historic flow of 1955 and 1790 are almost identical to that of the corresponding present field, whereas mean of NRM (after AF demagnetization) inclinations for each of the four pre-historic lava flow units is twelve to thirteen degrees lower than the present field inclination. We observed three cases of very large inclination variations from within a single flow, the best fitting curves of which are linear, second and third order polynomials each from within a single flow, whereas no present field variations are observed. This phenomena can be attributed to the notion that local magnetic anomalies on the surface of an active volcano are not permanent, but are transient. Therefore we believe that local magnetic anomalies of an active volcano may be constantly modified due to on going subsurface injections and circulations of hot material and also due to wide spacial and temporal distribution of highly magnetic basaltic flows that will constantly modify the topography which will in turn modify the local ambient geomagnetic field (Baag, et al., 1995). Our observations bring into question the general reliability of PSV data inferred from volcanic rocks, because on-going various geologic and geophysical activities associated with active volcano would continuously deflect and modify the ambient geomagnetic field.

• Journal of the Korean Geophysical Society
• /
• v.2 no.1
• /
• pp.75-88
• /
• 1999

A recent major subject of geophysical exploration is research into 3-D subsurface imaging with a composite information from the various geophysical data. In an attempt to interpret Schlumberger sounding data for the study area in 2-D and 3-D view, resistivity imaging was firstly performed and then pseudo-3-D resistivity volume was reconstructed by interpolating several 1-D resistivity plots. Electrical resistivity discontinuities such as fracture zone were successfully clarified in pseudo-3-D resistivity volume. The low resistivity zone mainly associated with fracture zone appears to develop down to granitic basement in the central part of the study area. Seismic velocity near the lineament is estimated to be approximately as small as 3,000 m/s, and weathering-layer for the southeastern part is interpreted to be deeper than for the northwestern part. Geophysical attributes such as electrical resistivity, seismic velocity, radioactivity for the Chojeong Area were analysed by utilizing a GIS software Arc/Info. The major fault boundaries and fracture zones were resolved through image enhancement of composite section (electrical resistivity and seismic refraction data) and were interpreted to develop in the southeastern part of the area, as characterized by low electrical resistivity and low seismic velocity. However, radioactivity attribute was found to be less sensitive to geological discontinuities, compared to resistivity and seismic velocity attributes.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.33-38
• /
• 2012

A resistivity method has been applied to wide range of engineering and environmental problems with the help of automatic and precise data acquisition. Thus, more accurate modeling and inversion of time-lapse monitoring data are required since resistivity monitoring has been introduced to quantitatively find out subsurface changes With respect to time. Here, we used the finite element method (FEM) for 3D resistivity modeling since the method is easy to realize complex topography and arbitrary shaped anomalous bodies. In the FEM, the linear elements, also referred to as first order elements, have certain advantages of simple formulation and narrow bandwidth of system equation. However, the linear elements show the poor accuracy and slow convergence of the solution with respect to the number of elements or nodes. To achieve the higher accuracy of finite element solution, high order elements are generally used. In this study, we developed a 3D resistivity modeling program using high order Serendipity elements. Comparing the Serendipity element solutions for a cube model with the linear element solutions, we assured that the Serendipity element solutions are more accurate than the linear element solutions in the 3D resistivity modeling.