• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.328-337
• /
• 2009

To investigate the relation between pore fluid conductivity and bulk resistivity of a rock sample it is assumed that electrolyte solution perfectly substitute the pore fluid that occupied the pore space within the sample in general. In this study, it is investigated that how much can the electrolyte solution substitute the pore fluid by repeating the same saturation process. Four kinds of NaCl solutions of 8, 160, 3200, 64000 ${\mu}S$/cm are used. The saturation process has repeated four times for each electrolyte in increasing conductivity order first then four times each in decreasing order. The more the saturation process repeated with the same electrolyte, the more electrolyte solution substitute the pore fluid. Geometric mean of bulk resistivity in increasing and decreasing orders with the same electrolyte solution is assumed to be mostly close to the bulk resistivity with perfect substitution. Bulk resistivity measurements for both increasing and decreasing order differs within 10% to the geometric mean when repeating the saturation process 4 times while maximum 40% difference is observed when single saturation process for each electrolyte solution with increasing order. The modified parallel resistant model can generally represent the relations between pore fluid resistivity and bulk resistivity in the experiment, but more experimental data with various rock samples with different porosity is needed to generalize the model.

• Geophysics and Geophysical Exploration
• /
• v.16 no.2
• /
• pp.71-78
• /
• 2013

The one-way wave equation migration is much more computationally efficient comparing with reverse time migration and it can provide better image than the migration algorithm based on the ray theory. We have developed the prestack depth migration module adopting (GS) propagator designed for vertical transverse isotropic media. Since GS propagator considers the higher-order term by expanding the Taylor series of the vertical slowness in the thin slab of the phase-screen propagator, the GS migration can offer more correct image for the complex subsurface with large lateral velocity variation or steep dip. To verify the validity of the developed GS migration module, we analyzed the accuracy with the order of the GS propagator for VTI media (GSVTI propagator) and confirmed that the accuracy of the wavefield propagation with the wide angles increases as the order of the GS propagator increases. Using the synthetic seismic data, we compared the migration results obtained from the isotropic GS migration module with the anisotropic GS migration module. The results show that the anisotropic GS migration provides better images and the improvement is more evident on steeply dipping structures and in a strongly anisotropic medium.

• Journal of the Korean Geophysical Society
• /
• v.8 no.1
• /
• pp.1-6
• /
• 2005

Defects of dam body which can be induced in seepage or leakage procedure can directly affect dam safety. Therefore, a proper inspection method should be carried out in the first place to find out their positions and sizes. After that, some reinforcement works such as grouting and the corresponding assessment could be taken in a proper way. The dam(center core type earth dam) issued in this study has been in need for intensive diagnosis and reinforcement work, because a lot of slumps similar to cracks, seepage and some boggy area have been observed on the downstream slope. High resolution seismic reflection method was performed on the crest profile twice before and after grouting work(Aug. 2001 and Nov. 2004) aimed at the dam inspection and the assessment of grouting efficiency as well. To enhance the data resolution, P-beam energy radiation technique which can reduce the surface waves and hence to reinforce the reflection events was used. Strong reflection events were recognized in the stack section before grouting work, It seems that the events would be caused by e.g. horizontal cracks with a considerable aperture. Meanwhile such strong reflection events were not observed in the section after grouting. That is, the grouting work was dear able to reinforce the defects of dam body. Hence, the section showed an well arranged picture of dam inner structure. In this sense, seismic reflection method will be a desirable technique for dam inspection and for monitoring dam inner structure as well.

• Journal of the Korean Geophysical Society
• /
• v.2 no.3
• /
• pp.191-200
• /
• 1999

Analysis of chirp high-resolution seismic profiles from the central Yellow Sea reveals that depositional environments in this area can be divided into three distinctive zones from west to east: (1) subaqueous delta system near the Shandong Peninsula, (2) erosional zone in the central Yellow Sea, and (3) tidal sand ridges and sand waves near the Korean Peninsula. The Shandong subaqueous delta, extending southward from the Shandong Peninsula, changes gradually into prodelta southeastward. The sediments originated from the Yellow River are transported southward along the Chinese coastal area. The erosional zone in the central Yellow Sea contains numerous paleochannels and shows linear erosional features trending northwest-southeast. The erosional zone would be dominated by non-depositional or erosional processes during the Holocene. Tidal sand ridges and sand waves are well developed along the western coast of Korea. The residual sands, which were originally fluvial sediments at the sea-level lowstand, are interpreted as the result of winnowing process during the sea-level rise. Modern sand ridges generally migrates in a northeast-southwestern direction, which coincide with dominant tidal current direction.

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

To estimate the shear-wave velocity (${\nu}_s$ beneath the OTVZ seismic station on Ngauruhoe volcano in New Zealand, we calculated receiver functions (RFs) using 127 teleseismic data ($Mw{\geq}5.5$) with high signal-to-noise ratios recorded during November 11, 2011 to September 11, 2013. The genetic inversion algorithms was applied to 21 RFs calculated by the iterative time-domain deconvolution method. In the 1-D ${\nu}_s$ model derived by the inversion, the Moho is observed at a 14 km depth, marked by a ${\nu}_s$ transition from 3.7 km/s to 4.7 km/s. The average ${\nu}_s$ of the overlying crust is 3.4 km/s, and the average ${\nu}_s$ of a greater than 9-km thick low-velocity layer (LVL) in the lower crust is 3.1 km/s. The LVL becomes thinner with increasing distance from the station. Another LVL thicker than 10 km with ${\nu}_s$ less than 4.3 km/s is found in the upper mantle. Those LVLs in the lower crust and the upper mantle and the relatively thin crust might be related to the magma activity caused by the subducting Pacific plate.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.323-329
• /
• 2010

Three deconvolution methods were applied to stacked seismic data obtained to investigate gas-hydrates in the Ulleung Basin, East Sea: (1) minimum-phase spiking deconvolution, (2) minimum-phase spiking deconvolution using an averaged wavelet from all traces, and (3) deterministic deconvolution using a wavelet with phases computed from well-logs. We analyzed the resolving property of these methods for lithological boundaries. The first deconvolution method increases temporal resolution but decreases lateral continuity. The second method shows, in an overall sense, similar results to the spiking deconvolution using a minimum phase wavelet for each trace; however, it results in a more consistent and continuous bottom-simulating reflector (BSR) and better resolved sub-BSR reflectors. The results from the third method reveal more detailed internal structures of debris-flow deposits and increased continuity of reflectors; in addition, the seafloor reflection and the BSR appear to have changed to a zero-phase waveform. These properties help more precisely estimate the distribution and reserves of gas hydrates in the exploration area by improving analysis of facies and amplitude of the BSR.

• 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.11 no.2
• /
• pp.116-126
• /
• 2008

In gas hydrate survey, seismic amplitude and frequency characteristics play a very important role in determining whether gas hydrate exists. According to the variation of source frequency and scatterer size, we study seismic amplitude characteristics using elastic modeling applied at staggered grids. Generally speaking, scattering occurs in proportion to the square of source frequency and the scatterer volume, which has an effect on seismic amplitude. The higher source frequency is, the more scattering occurs in gas hydrate bearing zone. Therefore, BSR is hardly observed in high frequencies. On the other side, amplitude blanking zone and BSR is clearly observed in lower frequencies although the resolution is poor as a whole. Seismic reflections traveling through free-gas layer below gas hydrate bearing zone decay so severely a high frequency component that a low frequency term is dominant. Amplitude anomaly of BSR result from high acoustic impedance contrast due to free-gas, which is a very crucial factor to estimate gas hydrate bearing zone. Seismic frequency analysis is carried out using wavelet transform method that frequency component could be decomposed with time variation. In application of wavelet transform to the seismic physical experiments data, we can observe that reflections traveling through air layer, which corresponds to the free-gas layer, decay a high frequency component.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.361-368
• /
• 2007

In this paper, we suggested the relationship between resistivity of coastal groundwater and NaCl equivalent salinity for the quantitative interpretation the results of surface/borehole resistivity and electromagnetic data. 38 groundwater samples having electrical conductivity higher than about 1,000 ${\mu}S/cm$ were analyzed to derive the empirical relationship between groundwater resistivity and NaCl equivalent salinity. We used Schlumberger chart GEN-8 to convert ion concentration from hydrochemical analysis to the equivalent NaCl salinity, and the portable meter to measure the in situ electrical conductivity of groundwater samples. From the hydrochemical analysis, relationship between the groundwater resistivity $(R_w)$ and equivalent NaCl salinity (Eq_NaCl) is expressed as Eq_NaCl=$5935.3551{\times}R_w^{-1.0993}$, and relationship between the groundwater electrical conductivity (EC) and total dissolved solids (TDS) is expressed as TDS=0.721*EC. We believe these relationships are very useful to assess the seawater intrusion in western and southern coastal area.