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

• The Journal of Engineering Geology
• /
• v.3 no.3
• /
• pp.231-239
• /
• 1993

This paper presents a result on the determination of velocity distribution by a tomographic inversion of crosshole seismic traveltimes in transversely isotropic(aniso tropic) media. The crosshole traveltimes used in this study are synthetic ones computed by ray tracing for some models having isotropic and transversely isotropic velocity distributions. The traveltimes are inverted by a general ART and ansotropic ART which considers the transversely isotropic effect during inver sion. The aniotropic ART gives accurate velodty distributions of transversely isotropic and isotropic models, while the isotropic ART determines accurate velocities only for the isotropic model but inaccurate for the transversely isotropic one. Therefore, the anisotropic ART may be used in case where no information is known on the isotropy or transverse isotropy of a survey area.

• 한국신재생에너지학회:학술대회논문집
• /
• 2006.06a
• /
• pp.382-385
• /
• 2006

한국지질자원연구원은 1997년부터 새로운 에너지 자원으로 활용 가능성을 포함하고 있는 가스 하이드레이트를 조사하기 위해 동해 일원에서 탄성파탐사를 실시하고 있다. 탄성파 반사 자료로부터 가스 하이드레이트 부존여부를 확인하는 방법은 해저면과 평행하면서 위상이 반대로 나타나는 고진폭 반사파 BSR (Bottom Simulating Reflection)과 BSR상부에서의 진폭감소, 하부에서 진폭증가와 구간속도 감소 등을 들 수 있다. 여기에서는 가스 하이드레이트 탐사자료에 대한 일반자료처리와 함께 BSR을 포함하고 있는 탄성파 반사자료에 대해 코드 병렬화된 PSPI를 이용하여 깊이영역 구조보정을 실시하였다. 고용량 탐사자료로 구성된 탄성파 반사자료에 깊이영역 구조보정을 적용하기 위해서는 고성능 컴퓨터와 병렬처리 기술이 필요하다. PSPI(Phase Shift Plus Interpolation)법은 적은 컴퓨터 계산량과 효율성 그리고 주파수 영역에서 구조적으로 병렬화가 용이한 특성을 지니고 있어 구조보정에 많이 이용되고 있다. 여기에서는 MPI(Message Passing Interface)-LAM을 이용하여 병렬코드화된 PSPI를 개발하고 인공합성모델과 동해 가스 하이드레이트 깊이영역 구조보정에 적응하였다.

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.50-57
• /
• 2011

For effective and accurate prediction of overpressure in the Efomeh field, located in the Niger delta basin of Nigeria, integrated seismic and borehole analyses were undertaken. Normal and abnormal pore pressure zones were delineated based on the principle of normal and deviation from normal velocity trends. The transition between the two trends signifies the top of overpressure. The overpressure tops were picked at regular intervals from seismic data using interval velocities obtained by applying Dix's approximation. The accuracy of the predicted overpressure zone was confirmed from the sonic velocity data of the Efomeh 01 well. The variation to the depth of overpressure between the predicted and observed values was less than 10mat the Efomeh 01 well location, with confidence of over 99 per cent. The depth map generated shows that the depth distribution to the top of the overpressure zone of the Efomeh field falls within the sub-sea depth range of 2655${\pm}$2m (2550 ms) to 3720${\pm}$2m (2900 ms). This depth conforms to thick marine shales using the Efomeh 01 composite log. The lower part of the Agbada Formation within the Efomeh field is overpressured and the depth of the top of the overpressure does not follow any time-stratigraphic boundary across the field. Prediction of the top of the overpressure zone within the Efomeh field for potential wells that will total depth beyond 2440m sub-sea is very important for safer drilling practice as well as the prevention of lost circulation.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.43-48
• /
• 1998

Recently seismic tomography has been widely used to visualize subsurface structure for resource explorations and construction site evaluation. We studied a way to include fresnel zone concept in the conventional ray-based traveltime tomography. The algorithm developed uses the same order of computing time as the conventional traveltime to mography but incorporates the rigorous wavepath concept of wave-equation tomography. Some experiments to synthetic and real data show reasonable results compared to conventional ray-based traveltime tomography.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.139-147
• /
• 2006

We had conducted a numerical modeling to investigate seismic properties of gas hydrate with field parameters acquired over the East sea in 1998. We used a 2-D staggered grid finite difference method to generate synthetic elastic seismograms for multi-channel seismic survey, OBC (Ocean Bottom Cable) survey and VCS (Vertical Cable Seismic) survey. The results of this study showed that the method using staggered grid yielded stable results and could be used to seismic imaging. We could find out the high amplitude anomaly and the phase reversal phenomenon of reflection wave at interface between the gas hydrate layer and free gas layer such a BSR (Bottom Simulating Reflector) which is the evidence for existence of gas hydrate in seismic reflection data. And we computed the reflection coefficients at the incident angles corresponding to offset distance with the synthetic seismograms. The reflection coefficients acquired from the numerical modeling were nearly consistent with the reflection coefficient computed by Shuey's equation.

• Geophysics and Geophysical Exploration
• /
• v.23 no.2
• /
• pp.97-114
• /
• 2020

Recently, many studies have been actively conducted on the application of machine learning in all branches of science and engineering. Studies applying machine learning are also rapidly increasing in all sectors of seismic exploration, including interpretation, processing, and acquisition. Among them, fault detection is a critical technology in seismic interpretation and also the most suitable area for applying machine learning. In this study, we introduced various machine learning techniques, described techniques suitable for fault detection, and discussed the reasons for their suitability. We collected papers published in renowned international journals and abstracts presented at international conferences, summarized the current status of the research by year and field, and intensively analyzed studies on fault detection using machine learning. Based on the type of input data and machine learning model, fault detection techniques were divided into seismic attribute-, image-, and raw data-based technologies; their pros and cons were also discussed.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.99-104
• /
• 2009

We propose a waveform inversion method for SH-wave data obtained in a shallow seismic refraction survey, to determine a 2D inhomogeneous S-wave profile of shallow soils. In this method, a 2.5D equation is used to simulate SH-wave propagation in 2D media. The equation is solved with the staggered grid finite-difference approximation to the 4th-order in space and 2nd-order in time, to compute a synthetic wave. The misfit, defined using differences between calculated and observed waveforms, is minimised with a hybrid heuristic search method. We parameterise a 2D subsurface structural model with blocks with different depth boundaries, and S-wave velocities in each block. Numerical experiments were conducted using synthetic SH-wave data with white noise for a model having a blind layer and irregular interfaces. We could reconstruct a structure including a blind layer with reasonable computation time from surface seismic refraction data.

• Geophysics and Geophysical Exploration
• /
• v.11 no.3
• /
• pp.242-249
• /
• 2008

Recently, AVO analysis has been widely used in oil exploration with seismic subsurface section as a direct indicator of the existence of the gas. In the case of the deep reservoirs like the gas reservoirs in the East-sea, it is often difficult to observe AVO responses in CMP gathers even though the bright spots are shown in the stacked section. Because the reservoir becomes more consolidated as its depth deepens, P-wave velocity does not decrease significantly when the pore fluid is replaced by the gas. Thus the difference in Poisson's ratio, which is a key factor for AVO response, between the reservoir and the layer above it does not increase significantly. In this study, we analyzed the effects of Poisson's ratio difference on AVO response with a variety of Poisson's ratios for the upper and lower layers. The results show that, as the difference in Poisson's ratio between the upper and lower layers decreases, the change in the reflection amplitude with incidence angle decreases and AVO responses become insignificant. To consider the limitation of AVO responses shown in the gas reservoir in East-sea, the velocity model was made by simulation Gorae V structure with seismic data and well logs. The results of comparing AVO responses observed from the synthetic data with theoretical AVO responses calculated by using material properties show that the amount of the change in reflection amplitude with increasing incident angle is very small when the difference in Poisson's ratio between the upper and lower layers is small. In addition, the characteristics of AVO responses were concealed by noise or amplitude distortion arisen during preprocessing. To overcome such limitations of AVO analysis of the data from deep reservoirs, we need to acquire precisely reflection amplltudes In data acquisition stage and use processing tools which preserve reflection amplitude in data processing stage.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.1-14
• /
• 2008

We present a method for interpreting seismic records with arrivals and waveforms having characteristics which could be generated by extremely inhomogeneous velocity structures, such as non-typical oceanic crust, decollement at subduction zones, and seamounts in oceanic regions, by comparing them with synthetic waveforms. Recent extensive refraction and wide-angle reflection surveys in oceanic regions have provided us with a huge number of high-resolution and high-quality seismic records containing characteristic arrivals and waveforms, besides first arrivals and major reflected phases such as PmP. Some characteristic waveforms, with significant later reflected phases or anomalous amplitude decay with offset distance, are difficult to interpret using only a conventional interpretation method such as the traveltime tomographic inversion method. We find the best process for investigating such characteristic phases is to use an interactive interpretation method to compare observed data with synthetic waveforms, and calculate raypaths and traveltimes. This approach enables us to construct a reasonable structural model that includes all of the major characteristics of the observed waveforms. We present results here with some actual observed examples that might be of great help in the interpretation of such problematic phases. Our approach to the analysis of waveform characteristics is endorsed as an innovative method for constructing high-resolution and high-quality crustal structure models, not only in oceanic regions, but also in the continental regions.