• 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.2 no.4
• /
• pp.184-190
• /
• 1999

A study of gas hydrate is a worldwide popular interesting subject as a potential energy source. A seismic survey for gas hydrate have performed over the East sea by the KIGAM since 1997. General indicators of natural submarine gas hydrates in seismic data is commonly inferred from the BSR (Bottom Simulating Reflection) that occurred parallel to the see floor, amplitude decrease at the top of the BSR, amplitude Blanking at the bottom of the BSR, decrease of the interval velocity, and the reflection phase reversal at the BSR. So the seismic data processing for detecting gas hydrates indicators is required the true amplitude recovery processing, a accurate velocity analysis and the AVO (Amplitude Variation with Offset) analysis. In this paper, we had processed the field data to detect the gas hydrate indicators, which had been acquired over the East sea in 1998. Applied processing modules are spherical divergence, band pass filtering, CDP sorting and accurate velocity analysis. The AVO analysis was excluded, since this field data had too short offset to apply the AVO analysis. The accurate velocity analysis was performed by XVA (X-window based Velocity Analysis). This is the method which calculate the velocity spectrum by iterative and interactive. With XVA, we could determine accurate stacking velocity. Geobit 2.9.5 developed by the KIGAM was used for processing data. Processing results say that the BSR occurred parallel to the sea floor were shown at $367\~477m$ depths (two way travel time about 1800 ms) from the sea floor through shot point 1650-1900, the interval velocity decrease around BSR and the reflection phase reversal corresponding to the reflection at the sea floor.

• 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.3 no.2
• /
• pp.53-60
• /
• 2000

The server/client systems using the web protocol and distribution computing environment by network was applied to the MT data processing based on the Java technology. Using this network based system, users can get consistent and stable results because the system has standard analysing methods and has been tested from many users through the internet. Users can check the MT data processing at any time and get results during exploration to reduce the exploration time and money. The pure/enterprised Java technology provides facilities to develop the network based MT data processing system. Web based socket communication and RMI technology are tested respectively to produce the effective and practical client application. Intrinsically, the interpretation of MT data performing the inversion and data process requires heavy computational ability. Therefore we adopt the MPI parallel processing technique to fit the desire of in situ users and expect the effectiveness for the control and upgrade of programing codes.

• 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.9 no.3
• /
• pp.225-230
• /
• 2006

Fluid replacement and conductivity logging have been applied to three boreholes in coastal aquifer in order to identify permeable fractures and to estimate the origin of saline groundwater. Fluid replacement technique measures and monitors the change of borehole fluid conductivity with depth under ambient or pumping condition after replacing the original borehole fluid with different one (by pumping out original one and injecting simultaneously new one at the hole bottom). After the replacement of borehole fluid, the change of fluid conductivity can be the direct indicator of the intake flow of formation water through aquifer such as permeable fractures or porous formations. The conductivity profiles measured with times therefore indicate the locations of permeable zone or fractures within the open hole or the fully slotted casing hole. As a result of fluid conductivity logging for three boreholes at coastal area in Yeonggwang, Jeonam Province, it is interpreted that the seawater intrusion in this area is not by remnant saline groundwater after land reclamation but mainly by intrusion of saline water through fractured rock. This approach might be useful for assessing the characteristics of seawater intrusion, the design of optimal pumping, the mitigation of seawater intrusion using freshwater injection, and estimating the hydraulic characteristics in coastal aquifer.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.255-262
• /
• 2009

Conventionally, vertical-seismic-profiling (VSP) survey that provides high-resolution information has mainly performed to obtain the exact depth of the gas hydrate-bearing sediment, which is one of the key factors in the development of the gas hydrate. In this study, we extracted interval velocities and created corridor stacks from the first domestic zero-offset VSP data, which were acquired with three component receivers at UBGH09 borehole in Ulleung Basin where gas hydrate exists. Then we compared the corridor stacks with a CMP stacked section from surface seismic data. First of all, we converted the signals recorded with three component receivers to true vertical and horizontal components by phase rotation, and divided the data into direct waves and reflected waves by wavefield separation processing. The trend of the interval velocity extracted from the zero-offset VSP was similar to that of the sonic log obtained at the same borehole. Because the interval velocity of the gas hydrate-bearing sediment above the BSR was high, and it decreased suddenly through the BSR, we could infer that free gas is accumulated below the BSR. The results of comparing the corridor stacks to the CMP stacked section of the surface seismic data showed that most reflection events agreed well with those in the surface CMP stacked section and that the phase-rotated VSP data corresponded better with the surface seismic data than the VSP data without phase rotation. In addition, by comparing a corridor stack produced from the transverse component with the CMP stacked section of the surface seismic data, we could identify PS mode-converted reflections in the CMP stacked section.

• Geophysics and Geophysical Exploration
• /
• v.12 no.3
• /
• pp.239-245
• /
• 2009

To verify the reinforcing effect of grouting materials composed of colloid cement and ordinary portland cement on the water leakage region in a small scale dike, we performed a tubecasing method and applied surface electrical resistivity survey including electrical resistivity tomography (ERT) to find resistivity variation before and after grouting. Hydraulic conductivities after grouting show 10 times lower than those of before grouting. These variation indicates that the cement grout blocks the leakage pathway effectively. As the results of dipole-dipole resistivity survey along the dike, resistivity distribution after grouting did not represent noticeable spatial variation in time. Resistivity monitoring results at the dike with vertical electrical sounding (VES) showed that the region of decreasing apparent resistivity was occupied by the grout after grouting. Predicted resistivities from the inversion of ERT data well matched with results of VES at the same regions. From the ERT using check holes to inspect the effect of grouting, we could find that the ERT is quite effective to identify spatially the grout region in a dike.

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