• Geophysics and Geophysical Exploration
• /
• v.21 no.1
• /
• pp.41-53
• /
• 2018

In the traditional seismic processing, multiple reflections are treated as noise and therefore they are eliminated during data processing. Recently, however, many studies have begun to consider multiples as signals rather than noise for seismic imaging. Multiple reflections can illuminate an area where primary reflections are not able to cover, thus it is allowed that a smaller number of shots and receivers are used for imaging large areas. In order to verify this, surface-related multiples were used for reverse-time migration (RTM), and then we compared the results with conventional RTM images which are generated from primary reflections. To utilize multiples, we separated multiples from whole seismic data using surface-related multiple elimination (SRME) method. Numerical examples confirmed that the migration using multiples can image wider area than the conventional migration, particularly in the shallow subsurface layers. In addition, the migration of multiples could eliminate the acquisition footprints.

• 한국해양학회지
• /
• v.21 no.2
• /
• pp.92-100
• /
• 1986

Laboratory investigations of physical (density and porosity), acoustic (velocity and velocity anisotropy), and electrical (resistivity and resistivity anisotropy) properties in deep-sea carbonate sequences at DSDPsites 288 and 289 in the western equatorial Pacific were made and correlated as a function of diagenesis. Profile of resistivity shows almost a mirror image of velocity indicating that electrical resistivity can be a useful variable to determine the diagenetic stage. Some fluctuations in acoustic and electrical properties near the zones of cherty and siliceous limestones for both sites imply significant changes in pore geometry due to interbedded silica. The significantly reduced pore throat size by the presence of silica which provides excess calcium carbonate to adjacent pore spaces is partly responsible for several jumps in acoustic and electrical properties of the zones. These observed geophysical data are interpreted as the result of silica diagenesis influencing carbonate diagenesis.

• The Journal of Engineering Geology
• /
• v.15 no.1
• /
• pp.19-27
• /
• 2005

The Sambo mine is located in Hae-je Myeon, Moo-an Gun, Chollanamdo, which consists of host gneiss and rhyolite possessing quartzite veins with other compositions such as gold, silver, and sublimated sulfur. The ore grade estimated from the core was 0.05～10.9g/t or less in gold and 0.05～389g/t or less in silver, indicating a partial mineralization. The purpose of this paper is to understand the subsurface structures and the distribution of mineralized bodies in the Sambo mine using a combined method of Schlumberger, Wenner, and Dipole-di-pole resistivity surveys on the surface and the resistivity tomography survey in boreholes. The result of three-dimensional resistivity inversion showed that the mineralized body is extended to 240m long in the N10°～20°E direction, with 30m wide and 80 m thick from the surface. The low resistivity zones (<1,000ohm-m) determined from the resistivity image were in good agreement with the mineralized bodies and weak zones identified from the logged cores.

• Geophysics and Geophysical Exploration
• /
• v.21 no.2
• /
• pp.94-102
• /
• 2018

To characterize virtual reflection images of deep subsurface by the method of seismic interferometry, we analyzed effects of offset range, ambient noise, missing data, and statics on interferograms. For the analyses, seismic energy was simulated to be generated by a 5 Hz point source at the surface. Vertical components of particle velocity were computed at 201 sensor locations at 100 m depths of 1 km intervals by the finite difference method. Each pair of synthetic seismic traces was cross-correlated to generate stacked reflection section by the conventional processing method. Wide-angle reflection problems in reflection interferometry can be minimized by setting a maximum offset range. Ambient noise, missing data, and statics turn to yield processing noise that spreads out from virtual sources due to stretch mutes during normal moveout corrections. The level of processing noise is most sensitive to amplitude and duration time of ambient noise in stacked sections but also affected by number of missing data and the amount of statics.

• 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.10 no.4
• /
• pp.259-267
• /
• 2007

An efficient algorithm for inverting static-shifted magnetotelluric (MT) data has been proposed to produce a three-dimensional (3D) resistivity model. In the Gauss-Newton approach, computational costs associated with construction of a full sensitivity matrix usually make 3D MT inversion impractical. This computational difficulty may be overcome by using approximate sensitivities. We use four kinds of sensitivities in particular orders in the inversion process. These sensitivities are computed 1) analytically for an initial, homogeneous earth, 2) exactly for a current model, 3) approximately by the Broyden method, and 4) approximately using the previous adjoint fields. Inversion experiments with static-shifted synthetic and field MT data indicate that inversion results are highly dependent on characteristics of data and thus applying various combinations of sensitivities is helpful in obtaining a good image of the subsurface structure with reasonable computation time.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.171-179
• /
• 2021

We outline a process for estimating Cole-Cole parameters from time-domain induced polarization (IP) data. The IP transients are all inverted to 2D Cole-Cole earth models that include resistivity, chargeability, relaxation time, and the frequency exponent. Our inversion algorithm consists of two stages. We first convert the measured voltage decay curves into time series of current-on time apparent resistivity to circumvent the negative chargeability problem. As a first step, a 4D inversion recovers the resistivity model at each time channel that increases monotonically with time. The desired intrinsic Cole-Cole parameters are then recovered by inverting the resistivity time series of each inversion block. In the second step, the Cole-Cole parameters can be estimated readily by setting the initial model close to the true value through a grid search method. Finally, through inversion procedures applied to synthetic data sets, we demonstrate that our algorithm can image the Cole-Cole earth models effectively.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.11 no.4
• /
• pp.361-371
• /
• 2009

Resistivity tomography is very effective geophysical method to find out the resistivity distribution and its change in time around a tunnel. Thus, the resistivity tomogram can provide helpful information which is necessary for the effective maintenance of the tunnel. However, an air filled tunnel severely distorts tomography data, especially when the current or potential electrode is placed near the tunnel. Moreover, the distortion can often lead to misinterpretation of tomography monitoring data. To solve these problem, we developed a resistivity modeling and time-lapse inversion program which include a tunnel. In this study, using the developed program we assured that the inversion including a tunnel gives much more accurate image around a tunnel, compared with the conventional tomogram where the tunnel is not included. We also confirmed that the time-lapse inversion of resistivity monitoring data defines well resistivity changed areas around a tunnel in time.

• Geophysics and Geophysical Exploration
• /
• v.26 no.4
• /
• pp.199-210
• /
• 2023

The physical properties of rocks in reservoirs change after CO₂ injection, we modeled a reservoir with a transition zone within which the physical properties change linearly. The function of the Wolf reflection coefficient consists of the velocity ratio of the upper and lower layers, the frequency, and the thickness of the transition zone. This function can be used to estimate the thickness of a reservoir or seafloor transition zone. In this study, we propose a method for predicting the thickness of the transition zone using deep learning. To apply deep learning, we modeled the thickness-dependent Wolf reflection coefficient on an artificial transition zone formation model consisting of sandstone reservoir and shale cap rock and generated time-frequency spectral images using the continuous wavelet transform. Although thickness estimation performed by comparing spectral images according to different thicknesses and a spectral image from a trace of the seismic stack did not always provide accurate thicknesses, it can be applied to field data by obtaining training data in various environments and thus improving its accuracy.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.302-309
• /
• 2008

ERT imaging, especially 3-D method, is a very powerful means to obtain a very high resolution image of the subsurface for geotechnical or hydrogeological problems. In this paper, we introduce two examples of successful case histories, where the imaging targets were three-dimensional. First example is the case of 3-D fracture imaging for hydrogeologic application. In this example, the borehole deviation was a critical problem in the ERT imaging and we could obtain real 3-D attitude of fracture system by including the borehole deviation in the inversion. In the second case, we did field experiment to image the empty tunnel with the size of $2m{\times}2m$ and the target was very clearly imaged in 3-D space. In these examples, we could show that 3-D ERT imaging is a very powerful tool for the 3-D subsurface imaging and the method can provide enhanced imaging capabilities especially for the 3-D targets such as fractures and cavities or tunnel.