Geobit을 이용한 가스 하이드레이트 탐사자료 처리

Seismic Data Processing For Gas Hydrate using Geobit

가스 하이드레이트는 전세계적으로 새로운 에너지 자원으로 활용 가능성을 포함하고 있어 연구가 활발하게 진행되고 있다. 한국자원연구소에서는 1997년 부터 동해에서 메탄 하이드레이트 부존 잠재력 규명을 위한 탄성파 탐사를 하고 있다. 탄성파 자료에서 하이드레이트 부존을 의미하는 일반적인 특성은 해저면과 평행하게 나타나는 BSR(Bottom Simulating Reflection)과 BSR 상부에서 보이는 진폭감소 그리고 BSR 하부에서 보이는 진폭증가와 구간속도의 감소 그리고 BSR에서 반사파의 역전현상 등이 있다. 따라서 위와 같은 하이드레이트 부존특성을 탐지하기 위한 목적으로 실시되는 자료처리는 진 진폭을 유지하는 자료처리, 정밀 속도분석 및 AVO분석 등이 요구된다. 본 연구는 1998년 동해에서 취득된 탄성파 탐사자료를 처리하여 하이드레이트 부존 가능성을 확인하고자 하였다. 적용된 자료처리 공정은 구형확산 보정과 주파수 필터링, 공심점 분류, 정밀 속도분석 공정 등이다. AVO분석은 이용된 현장자료가 AVO를 분석할 정도의 입사각을 유지하고 있지 않아 제외하였다. 정밀 속도분석은 반복적으로 속도 스펙트럼을 구하는 방법으로 정확한 중합속도 결정이 가능한 XYA를 이용하였으며 자료처리의 모든 공정은 국내고유의 탄성파 자료처리 소프트웨어인 Geobit 2.9.5 를 이용하였다. 자료처리 결과 음원위치 $1650\~1900$에서 해저면으로 부터 약 $367\~477m$ 깊이(왕복주시 약 1800ms)에 해저면과 평행하게 발달한 BSR을 확인할 수 있었으며, BSR부근에서 구간속도 감소 뿐만 아니라 해저면 반사파의 위상과 반대인 반사파 역전현상도 확인할 수 있었다.

A study of gas hydrate is a worldwide popular interesting subject as a potential energy source. A seismic survey for gas hydrate have performed over the East sea by the KIGAM since 1997. General indicators of natural submarine gas hydrates in seismic data is commonly inferred from the BSR (Bottom Simulating Reflection) that occurred parallel to the see floor, amplitude decrease at the top of the BSR, amplitude Blanking at the bottom of the BSR, decrease of the interval velocity, and the reflection phase reversal at the BSR. So the seismic data processing for detecting gas hydrates indicators is required the true amplitude recovery processing, a accurate velocity analysis and the AVO (Amplitude Variation with Offset) analysis. In this paper, we had processed the field data to detect the gas hydrate indicators, which had been acquired over the East sea in 1998. Applied processing modules are spherical divergence, band pass filtering, CDP sorting and accurate velocity analysis. The AVO analysis was excluded, since this field data had too short offset to apply the AVO analysis. The accurate velocity analysis was performed by XVA (X-window based Velocity Analysis). This is the method which calculate the velocity spectrum by iterative and interactive. With XVA, we could determine accurate stacking velocity. Geobit 2.9.5 developed by the KIGAM was used for processing data. Processing results say that the BSR occurred parallel to the sea floor were shown at $367\~477m$ depths (two way travel time about 1800 ms) from the sea floor through shot point 1650-1900, the interval velocity decrease around BSR and the reflection phase reversal corresponding to the reflection at the sea floor.