• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.80-83
• /
• 2001

Four seismic reflection horizons in 3-D seismic data, coherence derived from the seismic data, and 38 well logs from the Boonsville Gas Filed in Texas were tried to be integrated and visualized in 3 dimensions. Time surface was constructed from pick times of the reflection horizons. Average velocities to each horizon at 38 well locations were calculated based on depth markers from the well logs and time picks from the 3-D seismic data. The time surface was transformed to depth surface through velocity interpolation. Coherence was calculated on the 3-D seismic data by semblance method. Spatial distribution of the coherence is captured easily in 3-D visualization. Comparing to a time-slice of seismic data, distinctive stratigraphic features could be correctly recognized on the 3-D visualization.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.262-271
• /
• 2002

This study presents a practical procedure for the Bayesian inversion of geophysical data. We have applied geostatistical techniques for the acquisition of prior model information, then the Markov Chain Monte Carlo (MCMC) method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.291-298
• /
• 2002

We have developed a two-dimensional (2-D) magnetotelluric (MT) inversion algorithm, which can include topographic effects in inversion. We use the finite element method (FEM) to incorporate topography into forward calculation. Topography is implemented simply by moving nodes of rectangular elements in z-direction according to the elevation of air-earth interface. In the inversion process, we adopt a spatially variable Lagrangian multiplier algorithm in the smoothness-constrained least-squares inversion. The inversion algorithm developed in this study reconstructs subsurface resistivity structure quite well when topography variation exists. Also, it turns out to be effective in both resolution and stability from a model study and field data application.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.84-88
• /
• 2001

For the elastic migation, the velocity errors between the initial velocity model and true velocity model seriously affect the migrated images. The assumption of an initial velocity model, thus, is one of the critical factor for the successful migration. In case of applying the layered earth model as an initial velocity model, the layer boundary having large velocity contrast can not be defined well with conventional traveltime calculation algolithms and we have the difficulties for expressing the characteristics of the real subsurface. Smooth Background Model (SBM) we have applied as an initial velocity model in our study is characterized to be linearly varying the velocity with the depth, which can express the velocity variation in the subsurface properly. Thus it can properly be applied to traveltime calculation algolithms such as Vidale's method. In this study, Kirchhoff operator for prestack migration was used and the absolute amplitude obtained by modeling was applied as a weighted value to consider the true amplitude for initial model. Initial velocity model for migration was determined by using stacking velocity and we applied this model to real data.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.232-240
• /
• 2015

Magnetotelluric (MT) survey investigates electrical structure of subsurface by measuring natural electromagnetic fields on the earth surface. For the accurate interpretation of MT data, the precise three-dimensional (3-D) modeling algorithm is prerequisite. Since MT responses are affected by electrical anisotropy of medium, the modeling algorithm has to incorporate the electrical anisotropy especially when analyzing time-lapse MT data sets, for monitoring engineered geothermal system (EGS) reservoir, because changes in different-vintage MT-data sets are small. This study developed a MT modeling algorithm for the simulation MT responses in the presence of electrical anisotropy by improving a pre-existing staggered-grid finite-difference MT modeling algorithm. After verifying the developed algorithm, we analyzed the effect of vertical transversely isotropic (VTI) anisotropy on MT responses. In addition, we are planning to extend the applicability of the developed algorithm which can simulate not only the horizontal transversely isotropic (HTI) anisotropy, but also the tiled transversely isotropic (TTI) anisotropy.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.249-253
• /
• 2011

Based on intrinsic and scattering quality factor ($Q_i^{-1}$ and$Q_s^{-1}$) obtained from the seismic data of the southeastern Korean Peninsula, the expected coda quality factor (${Q_{Cexp}}^{-1}$) was theoretically calculated using multiple scattering model, and was compared with other quality factors such as $Q_i^{-1}$, $Q_s^{-1}$, and observed $Q_C^{-1}$ obtained by single scattering model. While the ${Q_{Cexp}}^{-1}$ values are typically comparable to the $Q_i^{-1}$ values, the $Q_C^{-1}$ values are different from the ${Q_{Cexp}}^{-1}$ values except for the higher frequency. Future works require to consider depth-dependent attenuation.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.23-32
• /
• 2012

The comparison between time-lapse seismic datasets is the most popular method in the reservoir monitoring. The method of extracting the changes only due to the change in the reservoir is the essential technique in the comparison of time-lapse seismic datasets. In the paper, the conventional cross-equalization approaches and an enhanced optimized approach have been tested and compared each other. As conventional approaches, the bandwidth equalization and phase rotation methods have been tested in frequency, time and mixed domains, respectively and their results were compared each other. In order to overcome the limit of the conventional approaches, which loses high frequency components, a new constrained optimum filtering method was proposed and experimented. The new constrained filtering method has shown the improvement in broadening the bandwidth of the components of reservoir changes by acquiring optimized match filter.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.214-219
• /
• 2011

We perform the frequency-domain waveform inversion based on the residual-selection strategy. In the residual-selection strategy, we classify time-domain residual wavefields into several groups according to the order of absolute amplitudes. Because the residual wavefields are normalized after regularization of the gradient directions within each group, the residual-selection strategy plays a role in enhancing the small-amplitude wavefields, which contributes to improving the deep parts of inverted subsurface images. After classifying residuals in the time domain, they are transformed to the frequency domain. Waveform inversion is performed in the frequency domain using the back-propagation technique which has been popularly used in reverse-time migration. The residual-selection strategy is applied to the SEG/EAGE salt and IFP Marmousi models. Numerical results show that the residual-selection strategy yields better results than the conventional frequency-domain waveform inversion.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.203-213
• /
• 2011

Time-lapse seismic surveys make repeated seismic surveys at different stages of oil production of a hydrocarbon reservoir to monitor changes in reservoir like fluid saturation. Since the repeatable surface seismic measurements can identify fluid types and map fluid saturations, oil and gas companies can make much more informed decision during not only production but also drilling and development. If time-lapse seismic surveys compare 3D seismic surveys, the time-lapse surveys are widely called as 4D seismic. A meaningful time-lapse interpretation is based on the repeatability of seismic surveys, which mainly depends on improved positioning and reduced noise (if surveys were designed properly through a feasibility study). The time-lapse interpretation can help oil and gas companies to maximize oil and gas recovery. This paper discusses about time-lapse seismic surveys mainly focused on feasibility, repeatability, data processing and interpretation.

• Geophysics and Geophysical Exploration
• /
• v.14 no.3
• /
• pp.191-202
• /
• 2011

Reserve-time migration (RTM) using a two-way wave equation is one of the most accurate migration techniques. RTM has been conducted by assuming that subsurface media are isotropic. However, anisotropic media are commonly encountered in reality. Conventional isotropic RTM may yield inaccurate results for anisotropic media. In this paper, we develop RTM algorithms for vertical transversely isotropic media (VTI) and tilted transversely isotropic media (TTI). For this, the pseudo-acoustic wave equations are used. The modeling algorithms are based on the high-order finite-difference method (FDM). The RTM algorithms are composed using the cross-correlation imaging condition or the imaging condition using virtual sources. By applying the developed RTM algorithms to the Hess VTI and BP TTI models, we could obtain better images than those obtained by the conventional isotropic RTM.