• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.45-54
• /
• 2001

A three-dimensional finite difference time-domain modeling algorithm based on staggered grid and considering transmitting and receiving antennas has been developed to simulate Ground Penetrating Radar (GPR) survey. This algorithm adopted the subcellular method to simulate the dipole antennas being used in GPR system and added resistors to reduce ringing caused by the reflections at the ends of an antenna. Comparison of the output voltages in the presence of the resistors for half-space said that the ringing and the amplitude of output voltage decreased as the number of resistors increased, and the antenna was designed based upon this result. Radiation patterns were derived to understand the distribution of electric field energy in the planes including or normal to the antenna. The electric field energy concentrated on vertical direction in the plane including antenna more than in normal plane.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.18-22
• /
• 2002

As a single arrival traveltime, maximum energy arrival traveltime has been known as the most proper operator for Kirchhoff migration. In case of the model having the simple structure, both the first arrival traveltime and the maximum energy arrival traveltime can be used as the correct operators for Kirchhoff migration. However for some model having the complex and high velocity contrast structure, the migration using the first arrival traveltime can't give the correct depth section. That is, traveltime to be required in Kirchhoff migration is the maximum energy traveltime, but, needs considerably more calculation time than that of first arrival. In this paper, we propose the method for calculating the traveltime approximated to the maximum energy arrival using one-way wave equation. After defining the WAS(Wrap Around Suppression) factor to be used for calculating the first arrival traveltime using one-way wave equation as the function of lateral grid interval and depth and considering the delay time of source wavelet. we calculate the traveltime approximated to the maximum energy arrival. to verify the validity of this traveltime, we applied this to the migraion for simple structure and complex structure and compared the depth section with that obtained by using the first arrival traveltime.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.33-45
• /
• 2002

It was necessary to devise new techniques to overcome and enhance the resolution limits of traveltime tomography. Waveform inversion has been one of the methods for giving very high resolution result. High resolution image could be acquired because waveform inversion used not only phase but amplitude. But waveform inversion was much time consuming Job because forward and backward modeling was needed at each iteration step. Velocity-stress method was used for effective modeling. Resolution limits of imaging methods such as travel time inversion, acoustic and elastic waveform inversion were investigated with numerical models. it was investigated that Resolution limit of waveform inversion was similar tn resolution limit of migration derived by Schuster. Horizontal resolution limit could be improved with increased coverage by adding VSP data in cross hole that had insufficient coverage. Also, waveform inversion was applied to realistic models to evaluate applicability and using initial guess of travel time tomograms to reduce non-linearity of waveform inversion showed that the better reconstructed image could be acquired.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.174-179
• /
• 1999

Shear wave velocity, one of major elastic constants in the dynamic design for civil structures, is conventionally measured from downhole, crosshole or sonic logging tests. SASW (Spectral Analysis of Surface Waves) method, which overcomes the disadvantage of the in-hole tests, can evaluate subsurface stiffness nondestructively and nonintrusively through measuring surface waves on surface. In this paper, principles of the SASW method are briefly described and the results of various field tests, conducted to investigate the applicability of the method, are summarized. The SASW method was successfully applied in evaluating the effects of dynamic compaction at Inchon international airport site, applied in evaluating the integrity of the lining and sidewall at a testing tunnel located in Mabukri, and applied in detecting thickness of a concrete retaining wall. The results of field tests and the nondestructive and economical characteristics of the method show the promising future of the SASW method in civil engineering projects.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.19-24
• /
• 2000

Seismic reflection data was obtained by using 28 and 100 Hz geophones at the ground subsidence sites in an old coal mine area. Frequency spectrum of the geophone analyzed with offset revealed that 1) In the near offset ($1\~10m$), the signals in the 100 Hz geophone data contains higher frequency components (up to 300 Hz) than that of the 28 Hz (<200 Hz), 2) In the intermediate offset ($11\~39m$), although the 28 Hz geophone data showed very similar frequency characteristics as the near offset data, the 100Hz geophone data seemed to be contaminated by noise at high frequency zone (>200 Hz). In the far offset ($\geq40\;m$), the signals in both the 28 and 100 Hz geophone data are attenuated to noise level at high frequency Bone more than 150 Hz.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.25-32
• /
• 2000

The geophysical data management system was developed to meet both the increasing demands of geophysical data in the practical application of civil engineering, underground water survey, and environmental problems and needs for digital archive and quality control of geophysical data. The system for a data manager is developed under Client/Server (C/S) environment. This manager system is characterized by a relational geophysical database system using MS SQL-server, standardization of geophysical data format, the development of C/S interface program for Windows environment, and the development of transfer program module for the searched data. The system developed for a general user under the internet environment is characterized by Web service (URL:http//geophy.kigam.re.kr) and the development of plug-in module to visualize geophysical image data.

• Geophysics and Geophysical Exploration
• /
• v.2 no.4
• /
• pp.191-201
• /
• 1999

The interpretation of resistivity data has, so far, mainly been made under the assumption that the earth is of relatively simple structure and then using one or two-dimensional inversion scheme. Since real earth structure and topography are fully three-dimensional and very complicated In nature, however, such assumptions often lead to misinterpretation of the earth structures. In such situations, three-dimensional inversion is probably the only way to get correct image of the earth. In this study, we have developed a three-dimensional inversion code using the finite element solution for the forward problem. The forward modeling algorithm simulates the real field situation with irregular topography. The inverse problem is solved iteratively using the least-squares method with smoothness constraint. Our inversion scheme employs ACB (Active Constraint Balancing) to enhance the resolving power of the inversion. Including Irregular surface topography in the inversion, we can accurately define the earth structures without artifact in the numerical tests. We could get reasonable image of earth structure by Inverting the real field data sets taken over highway bridge construction site.

• Geophysics and Geophysical Exploration
• /
• v.3 no.2
• /
• pp.70-75
• /
• 2000

In the seismic migration, Kirchhoff and reverse time migration are used in general. In the reverse time migration using wave equation, two-way and one-way wave equation are applied. The approach of one-way wave equation uses approximately computed downward continuation extrapolator, it need tess amounts of calculations and core memory in compared to that of two-way wave equation. In this paper, we applied one-way wave equation to pre-stack reverse time migration. In the frequency-space domain, forward propagation of source wavefield and back propagration of measured wavefield were executed by using monochromatic one-way wave equation, and zero-lag cross correlation of two wavefield resulted in the image of subsurface. We had implemented prestack migration on a massively parallel processors (MPP) CRAYT3E, and knew the algorithm studied here is efficiently applied to the prestck migration due to its suitability for parallelization.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.185-192
• /
• 2006

Cares should be taken when performing the P and S wave velocity loggings in engineering and environmental fields. Some of them are the effect of casing, which is installed to prevent the borehole collapsing when the drilling is done on the loose ground such as soil and/or soft rock, and the discrepancy of the velocities of the same media according to the difference of the source wave frequency spectrum. The elastic moduli obtained from the P and S wave velocity logging have the dynamic characteristics. To overcome these difficulties, the following suggestions are recommended; (1) develop and apply a careful drilling technique that can keep the borehole wall without a casing, and (2) apply the logging methods with the suitable frequency bandwidth for the object of the velocity logging. It is important to make the aseismological engineers understand the difference between the dynamic elastic moduli and the static ones obtained from mechanical test, and to advise them to use the information properly.